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SUMMARY 


Considerable  information  is  available  to  date  concerning  the  charact¬ 
eristics  of  the  flow  in  the  planetary  boundary  layer,  which  comprises  roughly 
the  lowest  two  thousand  feet  of  the  atmosphere.  Unfortunately,  in  many  in¬ 
stances,  the  results  frcm  different  sources  do  not  always  agree.  This  is  in 
part  due  to  the  extremely  complex  nature  of  the  flow  and  as  a  result,  some 
confusion  exists  as  to  the  exact  description:  of  the  planetary  layer.  For 
this  reason,  a  fairly  extensive  survey  of  the  existing  data  has  been  carried 
out,  and  the  results  of  this  review  are  presented  herein.  The  planetary  layer 
is  described  in  detail  with  respect  to  both  mean  velocity  and  turbulence,  and 
the  effect  of  thermal  stability  and  surface  conditions  oh  these  characteristics 
is  discussed.  Finally,  a  simplified  analytical  representation  of  the  flow  in 
the  planetary  layer  is  presented. 
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NOTATION 


A,  B,  C 


cu(n) 

E^Ck,  n) 

e 

f 

P 

f(p)»fu(px)»  etc* 


g 

g(r),  gy(rx) >  etc. 

i 

i 

1 

J 

K 


k 


n,  n' 


constants  relating  fluctuating  component  variances  to 
friction  velocity  (Sec.  3.5) 

constants  used  in  cross-correlation  spectrum  expressions 
(Sec.  3.6.5) 

co-spectral  density  function  (Sec.  2.^) 
four-dimensional  spectral  density  tensor  (Sec.  2.k) 
total  mean  turbulent  kinetic  energy 


=  k.z 


longitudinal  correlation  coefficients  (Sec.  2.1  and 
2.5) 


gravitational  acceleration 

lateral  correlation  coefficients  (Sec.  211  and  2.5) 
unit  vector  in  x-direction 

L  ;  also,  an  index  referring  to  u,  v,  or  w 
unit  vector  in  y- directi  on 
an  index  referring  to  u,  v,  or  w 

I  I 

von  K&mSn's  constant  **  0.4 
eddy  viscosity 

reduced  frequency  or  wave  number  =  n/U  cycles  per  foot 

value  of  k  at  maximum  value  of  the  non-dimensionalized 
spectrum  of  the  * i *  velocity  component 

wave  number  vector  s  i  k  +jk  +  kk;  also,  unit 
vector  in  z  -direction  X  ^ 

I  k| 

Prandtl  mixing  length 

length  scale  defined  by  Davenport  and  Harris  (Secs. 
3.6.5  and  3-7) 


0t  (du/dz)T 

scaling  length  defined  by  — 


integral  length  scales  of  turbulence  (Sec.  2.2  and  2.5) 
frequencies,  cycles  per  second 


vi 


p 

pressure 

«1J (n) 

quadrature  spectrum  (Sec.  2.4) 

r 

spatial  separation  vector  *ir  +  j  r  +kr;|r|sr 

*  x  y  -  z  - 

Ri 

gradient  Richardson  Number  (Sec. 3*2) 

Rij^rx,ry,rz,T^ 

general  double  velocity  correlation  tensor  (Sec.  2.1) 

Rij,Rij^rx^  etc* 

one-dimensional  correlation  tensors  (Sec.  2.1) 

*iJ<T) 

Eulerian  time-delay  correlation  tensor  (Sec.  2.1) 

t 

time 

T 

temperature 

V 

integral  time  scale  of  turbulence  (Sec.  2.2) 

U(t) 

total  velocity  component  in  x-direction  a  U  +  u(t) 

u 

mean  velocity  in  x-direction 

U(t) 

turbulence  velocity  vector  *  i.  u(t)  +  Jv(  t)  +  k  v(t) 

0 

T 

friction  velocity 

°G 

gradient  wind  velocity  (Sec.  3.1) 

u(t),  v(t),  w(t) 

fluctuating  velocity  components  in  x,  y  and  z  directions, 
respectively 

u'  ,  v' ,  w' 

rms  values  of  fluctuating  velocity  components; 
u'  =  JvF,  v's  w'  = 

X,  y,  Z 

Cartesian  coordinate  axes 

z 

height  above  earth's  surface 

z 

o 

roughness  length  (Sec.  3*4) 

ZG 

gradient  height  (Sec.  3*1) 

a 

mean  velocity  power  law  index  (Sec,  3*4);  also,  angle 
between  mean  wind  and  geobtrop'ic  direction 

7 

lapse  rate  =  -dT/dz 

yo 

constant  defined  in  Sec.  3*6.4 

.  (  n) 

coherence  function  (Sec.  2.4) 

’,ij(k’yl’)r2)>?ij(k’ry) 

magnitude  of  cross-correlation  spectrum  (Sec. 3. 6. 5) 

r 

dry  adiabatic  lapse  rat;  **  1°C/100  m  «  5.5°F/lOOO  ft. 

A 

wavelength  =  l/k  =  U/n,  feet 

vii 

n 


P 

0 


4lJ(n),0lJ(kx)’  etc- 

*o« 


T 

T 

O 

CO 

a 

a 

a 


air  density 

potential  temperature  s  temperature  which  a  volume  of 
air  assumes  when  brought  adiabatically  from  its  existing 
pressure  to  that  at  the  earth's  surface 

one-dimensional  spectral  density  functions  (Sec.  2.b) 

one-dimensional  spectral  density  functions  (Eq.  (29)) 

power  spectrum  of  total  velocity  component  in  x-cirection 

function  related  to  cross-correlation  spectrum  (Sec. 

3.6.5) 

incremental  time-delay;  also  shear  stress 

surface  shear  stress 

frequency,  radians  per  second,  ■  2m 

reduced  frequency,  radians  per  foot,  *  2 77k 

non-dimensional  form  of  a 

time  average  of  a 
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I .  INTfoUKTTON 


In  the  past  few  years,  it  has  been  increasingly  important  to  know  in 
detail  the  characteristics  of  the  atmosphere  at  low  altitudes.  The  present  empha¬ 
sis  on  low-flying  and  V/STOL  aircraft  is  in  part  responsible  for  this  require¬ 
ment,  as  is  the  advent  of  rocketry,  with  its  large  launch  vehicles.  In  addition, 
the  continuation  of  the  trend  toward  taller  and  more  radically  shaped  buildings 
and  structures  necessitates  a  good  simulation  of  the  flows  at  these  levels  in  order 
realistically  to  determine  their  response  to  these  flows.  The  prominence  of  the 
pollution  control  issue  in  recent  years  has  also  added  weight  to  the  requirement 
of  achieving  a  proper  description,  and  ultimately  simulation  of  the  planetary 
boundary  layer. 

To  date,  considerable  data  has  been  gathered  in  both  the  aeronautical 
and  meteorological  disciplines.  While  the  information  is  ,by  no  means  complete, 
it  is  possible  to  draw  many  conclusions  about  the  nature  of  the  planetary  layer, 
and  in  particular  the  turbulence  therein.  Due  to  the  extremely  complex  nature 
of  the  flows,  however,  there  is  sometimes  considerable  variability  in  the  con¬ 
clusions  that  have  been  drawn.  For  example,  there  have  been  widely  scattered 
conclusions  about  the  magnitude  and  variation  of  the  integral  scales  of  turbu¬ 
lence  through  the  planetary  layer.  This  has  been  due  in  part  to  a  lack  of  rigor 
in  determining  exactly  what  scales  are  being  obtained  in  some  cases,  and  in  part 
to  the  refusal  of  the  atmosphere  to  submit  easily  to  simple  descriptions. 

It  is  the  purpose  of  this  review  to  describe  and  summarize  most  of 
the  existing  data  on  the  characteristics  of  the  flows  in  the  planetary  layer,  and 
to  attempt  to  clarify  some  of  the  discrepancies  that  have  occurred  in  the  past. 

In  addition,  in  light  of  the  data  available,  a  simplified  mathematical  model  is 
suggested  to  represent  the  flows. 


II.  DEFINITIONS  AND  THEORY  OF  TURBULENCE 


In  the  literature  on  turbulence  there  is  on  the  whole  a  fairly  wide 
range  of  notation  used.  In  order  to  avoid  confusion  and  to  clarify  precisely 
what  is  be^ng  referred  to  in  this  review,  the  notation  and  definitions  of  turbu¬ 
lence  theory  used  in  this  report  are  presented  in  some  detail  in  the  following 


pages . 
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Consider  the  two  arbitrary  points  P  and  P  separated  by  a  vector 
r  i  rv  ■*  J  rv  f  k  *n  a  homogenous  three-dimensional  turbulence  field  as  shown 

above.  Let  the  flow  be  moving  with  mean  velocity  0  along  the  x-axis  of  a  fixed 
Eulerian  reference  frame  of  which  PQ  is  the  origin.  The  points  P^,  P2  and  P^  are 
the  projections  of  P  on  the  coordinate  axes.  The  turbulent  motion  at  any  point 
in  the  field  is  represented  by  the  vector  U(t)  =  iu(t)  +  jv(t)  +  kw(t),  with 
appropriate  subscripts,  where  u(t),  v(t)  and  w(t)"are  the  longitudinal,  lateral 
and  vertical  fluctuating  velocity  components,  respectively.  Thus  the  total 
(Eulerian)  velocity  in  the  x-direction  at  PQ  is 

UQ(t)  =  0(t1,  T)  +  uQ(t) 

where  U( t.^.,  T)  is  the  mean  velocity  and  is  defined  by 

tx+  T 

U(ti,  T)=^  J  UQ(t)  dt  (1) 

VT 

and  is  the  same  at  all  points  in  the  field  due  to  the  assumption  of  homogeneity. 
The  mean  of  the  fluctuating  component  uft(t)  as  defined  by  Eq.  (l)  is  necessarily 
zero,  as  it  is  for  the  other  two  components  vQ(t)  and  wQ(t)  at  Pq  since  there  is 

mean  motion  only  in  the  x-direction,  and  similarly  for  any  other  points  in  the 
field.  The  three  fluctuating  components,  in  addition  to  U  (t),  are  random  func¬ 
tions  of  time,  and  as  such  are  said  to  be  stationary  if  their  statistical  prop¬ 
erties  are  independent  of  time  in  the  limit  as  T  -»°o.  Thus  for  stationary  flow, 
U(t  ,  T)  is  independent  of  t^  and  the  mean  velocity  is  given  by 

U  =  lim 


and  thus 

UQ(t)  =  U  +  UQ(t). 

Stationarity  of  the  velocity  components  will  be  assumed  throughout  this  review. 

Any  pair  of  fluctuating  components  at  either  P  or  P  which  are  separa¬ 
ted  in  time  or  space  can  be  multiplied  and  averaged  to  form  a  (double)  velocity 
correlation  function.  The  complete  set  of  possible  correlations  is  given  by 


R.  .(r,  t) 
ij  -* 


Ri.1(rx>WT)  i»J  =  U’V’W 

U^t)  U  (t  +  t) 


(3) 


which  is  a  second  order  tensor  whose  components  are  functions  of  the  four 
variables  shown  and  which  is  usually  called  the  'general  correlation  tensor'.. 
As  an  example  of  one  of  the  components  of  this  tensor,  consider  i  =  u,  and 
j  =■  v.  Then 


2 


1 


[ 


R 

uv 


3  UQ(t)  V(t  +  T) 

i  rT 

=  Urn  ^  u  (t)  v(t  +  t)  dt 
T-»w> 


=  lim 
T-»  oo 


1_ 

2T 


rT 

I  u(0,0,0,t)  v(rx,ry,rz,t+  -r)dt 

-T 


since 


uQ(t)  is  the  component  at  Pq  and  v(t)  is  at  P. 


(M 


It  should  be  noted  that  in  defining  the  mean  velocity  and  correlations 
above,  time  averages  of  the  velocity  signals  have  been  used,  since  it  is  these 
averages  that  are  actually  measured  in  practice.  However,  it  may  in  some  cases 
be  mathematically  advantageous  to  consider  ensemble  averages  rather  than  time 
averages.  For  example,  if  IL(t)  represents  an  ensemble  of  N  records  of  the  velo¬ 
city  U(t)  where  K  =  1,...N,  £hen  the  ensemble  average  of  U(t)  is  defined  by 


U 


=  lim 
N  ~>oo 


N 

I  VV 

K=1 


where  U^t)  is  assumed  stationary  such  that  U  is  independent  of  the  arbitrary  point 
in  time,  t^.  In  the  special  case  when  the  time-averaged  mean  velocity  as  defined 
by  Eq.  (2),  as  well  as  any  other  statistical  properties  of  U(t),  are  identical 
regardless  of  which  of  the  K  records  of  U(t)  is  used  in  obtaining  them,  U(t)  is  said 
to  be  ergodic.  In  this  case,  time  averages  are  equivalent  to  ensemble  averages 
such  that  the  mean  velocity  as  defined  by  Eq.  (2)  is  identical  to  the  ensemble 
average  defined  above.  In  practice,  ergodicity  of  the  velocities  is  usually  assumed. 
Note  also  that  ergodicity  requires  stationarity,  but  not  vice-versa. 


Of  the  general  correlation  tensor  of  Eq.  (3),  we  are  particularly  int¬ 
erested  in  five  special  cases.  These  are 


,  R^j(0»ry,0,0) ,  R^j (0,0,rz,0) , 

R^j(0,0,0,0)  and  R^j(rx,0,0,  t) . 


The  first  three  of  these  special  tensors  are  zero  time-delay  spatial  cross-corre¬ 
lations  between  the  turbulence  components  at  P  and  those  at  P  ,  P?  and  P_, 
respectively,  as  shown  in  the  above  sketch.  T$ey  are  given  typically  by 


R . . (r  ,0,0,0)  = 
l  j  x 


and  similarly, 


and 


uQ(t)ul(t) 

U0(t)Vl(t) 

uQ(t)w1(t) 

v0(t)ui(t) 

v0(t)vx(t) 

vQ(t)wl(t) 

w0(t)ui(t) 

wo(t)Vi(t) 

RiJ(0>ry 

V0-0’ 

,0,0)  =  RtJ( 

v°> 5  V 

V 

rz>- 

(5) 


(6) 
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The  fourt.h  tensor  is  the  well-known  Reynolds  Stress  tensor  and  is 


given  by 


R.  . (0,0, 0,0)  = 


Uo^(t)  uo(t)vo(t) 

v0(t)u0(t)  vo^(t> 

w  (t)u  (t)  w  (t)v  (t) 
o  O  0  0 


u0(t)wo(t) 

voh)»0(t) 

v'w 


(7) 


The  diagonal  components  of  this  tensor  are  of  course  the  familiar  mean-square  values 
ity  cfljimonents .  The  root-mean- square  (rms)  values  or  variances 
and  /w0-  are  rewritten  as  uQ' ,  vq'  and  w^',  respectively,  for  simplicity. 

Finally,  the  fifth  special  tensor  is  important  because  it  in  fact  represents  the 
single  point  time-delay  correlations  that  are  measured  by  a  probe  at  a  fixed  point 
such  at.  F Q  in  a  laboratory  ( Euler ian)  frame.  That  is,  in  a  time  increment  t,  the 
mean  flow  moves  along  the  x-axis  a  distance  Ut  and  if  we  simply  set  r  =  Ut,  we 

obtain  x 


R 

wt.ere 


and  is  obtained  by  correlating  time-delayed  signals  as  measured  by  a  probe  fixed 
at  Pq.  The  diagonal  components  of  this  tensor  are  called  'autocorrelations'  since 
the  components  are  being  correlated  with  themselves.  It  is  also  to  be  noted  that 
from  Eq^8)  for  t  =  0,  <R.  .  (t)  may  be  related  to  the  Reynolds  Stress  tensor  of 

Eq,(7);  that  istf?„(0)  =  R^ . 


(r  ,0,0,t)  = 

1  A 

Rij(frr,0,0,T)  = 

«U<T> 

(8) 

u0(t)uo(t+T) 

uo(t)vo(t+T) 

Uo(t)wo(t+T) 

(9) 

vo(t)uo(t+T) 

Vo(t)vo(t+T) 

vo(t)wo(t+T) 

wo(t)Uo(t+T) 

w0(t)vo(t+T) 

W  (t)w  (t+T) 

0  0 

Generally,  the  various  correlation  functions  defined  above  are  normalized 
to  yield  correlation  coefficients.  This  normalization  is  done  using  the  appropriate 
rms  values  of  the  velocity  components  such  that,  for  example,  the  normalized  form 
of  the  u-v  correlation  function  defined  in  Ec^(4)  is 


uv 


R  (r  ,r  ,r  ,  t) 
/  \  _  uv  x  y’  z 

(r  ,r  ,r  ,t)  = - : — * - 


and  similarly  for  any  of  the  other  correlations.  The  only  exceptions  to  this 
~ule  are  the  three  diagonal  components  of  the  Reynolds  Stress  tensor.  Since 
these  components  are  themselves  the  mean-square  values  of  the  velocities, 
normalisation  as  above  necessarily  yields  unity  and  is  therefore  of  no  use.  For 
this  reason,  the  roots  of_these  components  are  simply_non-dimensionalized  either 
by  the  friction  velocity  U^,  or  by  the  mean  velocity  U  in  which  case  the  familiar 
turbulence  intensities  given  by  uq'/ U,  vq'/U  and  wq'/U  are  obtained. 


In  general,  the  diagonal  components  of  the  spatial  correlation  co¬ 
efficient  tensors  are  referred  to  as  either  'longitudinal'  or  'lateral'  corre¬ 
lation  coefficients.  This  nomenclature  derives  from  the  direction  of  the  velo¬ 
city  component  being  correlated  at  the  two  points  with  respect  to  the  vector  r 
. eparating  the  points.  If  the  component  is  parallel  to  this  vector,  the 
correlation  is  called  longitudinal  and  the  function  f(r)  is  used  to  denote  the 


4 


coefficiert,  where 

r  = 

if  it  is  perpendicular  to  r,  the  correlation  is  called  lateral  and  the  function  g(r) 
denotes  the  coefficient.  Thus  there  are  three  longitudinal  correlation  coefficients 


giver  by 


and 


f  (r  ) 
u'  x' 

f  (r  ) 

v  y 

f  (r  ) 
wv  z7 


Euu(rx>* 

Ew(ry>’ 

Ewv(rz>  > 


(10) 


a.nd  six  lateral  coefficients  given  typically  by 

g  (r  )  =  R  (r  ),  g  (r  )  =  R  (r  ),  (ll) 

°vx  x7  wx  x7’  x7  wwx  x7*  '  7 

and  similarly  for  g  (r  ) ,  g  (r  ) ,g  (r  )  and  g  (r  ). 

°uv  Y  w  y  u'  z7  v  z7 

2.2  Integral  Scales 


An  integral  or  ’macro-scale’  of  turbulence  can  be  defined  for  any 
correlation  coefficient  1  of  any  of  the  special  tensors  defined  in  the  previous 
section  except  the  Reynolds  Stress  tensor.  The  scale  is  defined  as  the  integral 
of  the  correlation  coefficient  over  the  positive  range  of  its  independent  vari¬ 
able,  and  it  can  thus  be  either  a  length  scale  or  a  time  scale.  Ojje  can  there¬ 
fore  define  four  scale  tensors  corresponding  to  the  four  tensors  R.  (r  ), 

R^(r  ),  Rij(rz)  and^j(x).  However,  only  the  diagonal  components1of  #hese  tensors 

are  o?  major  interest,  and  in  the  case  of  the  first  three  tensors,  one  obtains 
nine  scale  lengths  defined  by 

rs  00 


(rx)drx’ 


(a) 


and 


Liy  s/o  Wdy 

(b) 

(12) 

z  r 00 

L1  *  Kil(rz)dV  1  = 

J  6 

(c) 

From  the  f?.  ,(t)  tensor,  the  three  time  scales  of  importance  are  given  by 


TiX  =  ^ii(T)  dT  (13) 

J  o 

x 

and  these  time  scales  can  be  related  to  the  length  scales  Li  through  the  use  of 
Taylor's  Hypothesis  (see  below). 


As  with  the  correlation  coefficients  of  Eqs^lO)  and  (ll),  the  length 
scales  of  Eqs.(l2)  are  referred  to  as  longitudinal  or  lateral  scales,  the  criterion 
logically  being  from  which  type  of  correlation  coefficient  they  are  obtained. 

Thus  from  these  equations  it  is  seen  that  L  x,  L  ^  and  L  z  are  longitudinal 
scales,  while  the  other  six  scales  of  Eqs,(l^;  )  are  lateral  scales. 


2 . 3  ^aylor's  Hypothesis 

Taylor's  Hypothesis  provides  a  time- space  transformation  which  allows 
spatial  variations  in  the  turbulence  field  to  be  expressed  in  terms  of  time 
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variations  at  a  fixed  point  in  the  field.  The  hypothesis  states  that  in  a  Lagrangian 

reference  frame  moving  with  the  mean  flow,  there  is  no  time  variation  of  the  comp¬ 
onents  -  that  is,  the  field  is  'frozen'.  Mathematically,  this  may  be  stated  as 

S/St  =  -OS/Sx  (l4) 

and  the  hypothesis  is  made  under  the  assumption  that  0  »  u'  .  The  physical  inter¬ 
pretation  of  this  is  that  the  time  fluctuations  at  a  fixed  point  in  the  field  such 
as  P0  can  be  imagined  to  be  caused  solely  by  the  entire  field  being  frozen  at  a 
particular  instant  and  convected  past  the  point  with  the  constant  velocity  0.  The 
velocity  fluctuations  over  a  period  of  time  at  the  point  will  then  be  identical 
with  the  instantaneous  distribution  of  the  velocity  u  (t)  along  the  x-axis  through 
the  point.  These  physical  implications  will  be  further  discussed  in  the  following 
section . 


In  terms  of  the  previously  defined  correlations,  Taylor's  Hypothesis  means 
that  Rij(£>T)  is  now  replaced  by  R^(r).  Thus  the  time-delay  correlation  of  Eq(8) 

may  now  be  written 


-  RiJ<V0'0> 


which  is  identical  to  the  spatial  cross-correlation  R. .(r  )  of  Eqi5). 


Thus  we 


or 


*u(t)  'W-  rx  =  °T 


(15) 


and  R^j(rx)  may  be  determined  simply  by  measuring  the  velocity  signals  at  a  fixed 

point  and  obtaining  the  time-delay  correlations  of  these  signals.  Using  this  re¬ 
sult  and  Eqs,(l2a)  and  ( 13 ) ,  it  is  now  possible  to  relate  the  scales  L  x  and  T.  . 
That  is,  ll 


r 00  ~ 

=  /  R. . (r  )dr 
ll  x  x 

J  o 


<Ru(t)  d(Ur)  =  to* 


(16) 


Strictly  speaking,  Taylor's  Hypothesis  was  made  for  homogeneous,  iso¬ 
tropic  turbulence,  although  it  has  been  found  quite  reasonable  for  non-isotropic 
uniform  flows.  Its  application  in  shear  flows,  however,  requires  further  con¬ 
sideration  (see  Sec.  3.3*^) • 


2.1  Spectral  Density  Functions 


Associated  with  any  correlation  function  is  a  corresponding  spectrum 
function  defined  as  the  Fourier  Transform  of  the  correlation.  Thus  for  the 
general  correlation  tensor  of  Eq,(3),  one  obtains  a  tensor  Of  spectra  given  by 

<»  -i27r(k.r  +  n'T) 

EiJ(k,n’)  *  16  fJJJ  R^r.-Oe  dr  dr  (17) 

—  00 


where  _k  -  i_  kx  +  k  =  k  is  the  wave  number  vector  corresponding  to  the  sep¬ 
aration  vector  r  after  the  transformation.  Writing  Eq(l7)  out  in  full  gives 


E.  . ( k  ,k  ,k  ,n' )  =  16  I  I  R.  .(r  ,r  ,r  ,  T)e 
ij  x’  y5  z*  IJJ I  iy  x’  y*  z’ 


-i27r(k  r+kr  +  kr  +  n'T) 
xx  y  y  z  z 


dr  dr  dr  dT. 
x  y  z 

(18) 
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The  tensor  E^,(k,n')  is  a  tensor  of  four-dimensional  spectrum  functions  and  can  of 
course  be  used  to  regain  the  correlation  tensor  through  the  inverse  Fourier  Trans¬ 
form  such  that 


Vi’T) 


1 

15 


ffffE iAn,) 


i27r(k.r  +  n'r) 


dkdn' 


which  was 


Suppose  now  that  we  insert  in  the  above  equations  the  relation  r  =  Ut 
used  to  relate  r^  and  t  in  the  previous  section.  Then  Ec^(l8)  becomes 


E.  .[(k  h  '/U),k  ,k  ]  =  8 
ij1-'  x  '  y  zJ 


R,,(r  ,r  ,r  )e 
ij  x’  y’  z' 


-i27r[(k  +n'/U)r  +k  r  +k  r  ] 
L  x  '  *  x  yy  z  z J 


dr  dr  dr  . 
x  y  z 


Now  let  n  -  k^U  +  n'  .  Physically,  the  frequency  n  represents  the  overall  time 

fluctuations  as  seen  by  an  observer  at  a  fixed  point  such  as  P  while  the  flow 
sweeps  by  with  mean  velocity  U  relative  to  the  point.  It  is  therefore  the 
frequency  which  would  be  measured  in  an  experiment  by  a  probe  fixed  at  P  ,  and 

is  the  sum  of  the  two  frequencies  n'  and  k  U.  The  frequency  n'  is  that  which 
would  be  seen  by  the  observer  at  a  fixed  point  in  a  Lagrangian  frame  moving 
with  the  fluid.  In  this  case,  the  relative  velocity  between  the  observer  and 
the  flow  would  be  zero  such  that  n  =  n' .  The  remaining  frequency  k  U  represents 

the  time  fluctuations  that  are  due  entirely  to  the  flow  field  being  convected 

past  P0  with  the  relative  mean  velocity  U.  In  fact,  if  U  is  large  enough  such 
that  k  U  »n' ,  then  n  «  k  U  and  this  is  none  other  than  Taylor's  Hypothesis. 

That  is,  if  the  flow  were  truly  'frozen',  n'  would  be  zero  and  n  =  k  0.  Thus  it 

can  be  seen  that  the  validity  of  Taylor's  Hypothesis  depends  on  the  relative 
velocity  between  the  observer  and  the  flow  field,  and  if  this  velocity  is  very 
small  such  that  n'  is  of  the  order  of  kxU,  then  n'  cannot  be  neglected  with 

respect  to  k  U  and  the  frozen  flow  hypothesis  is  not  valid.  Skelton  (Ref. 6) 
suggests  tha£  the  minimum  relative  velocity  for  the  hypothesis  to  be  valid  is  about 
one-third  of  the  mean  flow  velocity. 


Returning  now  to  the  spectrum  tensor,  we  may  write 

«  -i27r[(n/U)r  +k  r  +k  r  ] 

rrr  x  y  y  z  zJ 

E.  ,{n/U.k  ,k  )  =8  / / /  R.  . (r  ,r  ,r  )e  dr  dr  dr 

'  3  y  z  JJJ  iJ  x’  y’  z'  x  y  z 

—  00 

and  if  Taylor's  Hypothesis  is  assumed  (n*  =0),  n  =  k  U  and  E.  .  may  be  written 

X  1 J 


E.  .(k  ,k  ,k  )=  8 

ijv  x’  y’  z ' 


R.  .(r  ,r  ,r  )e 
x’  y*  z 1 


-i27r(k  r  +  k  r  +  k  r  ) 
xx  y  y  z  z' 


dr  dr  dr  . 
x  y  z 


(19) 


Consider  now  the  special  tensors  defined  by  Eqs.(5),  (6)  and  (9).  By 
Fourier  transforming  these  one-dimensional  correlations,  one  can  define  the 
corresponding  one-dimensional  spectrum  functions  by 


p  00 

<t>.  .(k  )  s  2  /  R.  .(r  ) 

iJ  x'  J  iy  x' 


-i27rk  r 

XX, 

e  dr 


(a) 
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00 


♦nOO  s  2  R^(0 


-i27rk  r 

y  y 


y  y 


0.  .(k  )  = 

IJ  Z 


-oo 

oo 

—  oo 


ij  y 


dr 


R.  .(r  )  e 
ij  z 


■i27rtc  r 

z  z 


dr 


(b) 

(c) 


(20) 


and 


poo  -127mT 

0.  (n)  =  2  /  ^  (t)  e  dT. 


(21) 


The  corresponding  inverse  transforms  are  of  course 


p  00 

i2vk  r 

=s  1 
~  2 

j  hj(kx)e 
J  —00  u 

X  X 

dk 

X 

(a) 

p  oo 

i27rk  r 

=  l 
—  2 

/  0.  .(k  )e 

y  y 

dk 

(b) 

J-co  ^  y 

y 

r  00 

i27rk  r 

=  1 
~  2 

J  0iJ(kz)e 

U  — OO 

z  z 

dk 

z 

(c) 

(22) 


and 


V7)  a  kf 


hj(n)e 


i27nrr 


dn. 


(23) 


The  0ij(n)  spectra  of  Eq.(2l)  are  those  which  are  obtained  by  placing  a  measuring 
instrument  at  a  fixed  point  in  the  flow,  and  if  Taylor's  Hypothesis  is  assumed, 
these  can  be  easily  related  to  the  0ij(kx)  spectra.  That  is,  from  Eqs.(l5)  and 

(20a),  and  since  r  =  Ut  and  k  =  n/U, 


0 


p  W 

iA>  =2j  _  Ru(r*) 


=  211 


-i27fk  r 

x  x  , 

j  v *  _  /€  dr 

-i27mT 

(t)  e  dT. 


Comparing  this  with  Eq.  (2l),  it  can  therefore  be  seen  that 

V11*5  =“'Vn)- 

Notice  also  that  in  the  particular  case  of  t  -  0,  Eq.(23)  may  be  used  to  relate 
the  Reynolds  Stress  tensor  to  the  0^(n)  tensor  by 

»  00 

(24) 

1  .1 

—  00 


p  00 

V0*  ’  Rij  =  5  J  *U(")dn- 


The  one-dimensional  spectra  defined  by  Eqs.(2C)  could  also  have  been 

obtained  by  integration  of  the  three-dimensional  spectra  of  Eq.(l9).  For  example, 

0.  (k  )  is  aiso  given  by  the  double  integration  of  E. .(k  ,k  ,k  )  over  k  and  k  , 
ij  x  ij  x  y  z  y  z 

Similarly,  two-dimensional  spectra  can  be  defined  by  integrating  E.  .  over  only  one 
variable,  but  these  will  not  be  dealt  with  here.  J 

As  with  the  correlation  functions  and  integral  scales  of  turbulence, 
special  terminology  is  generally  applied  to  the  spectrum  functions.  Of  the  nine 
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diagonal  components  of  the  tensors  defined,  by  Eqs.  (20),  uhree  are  termed  longi¬ 
tudinal  spectra  and  the  other  six  lateral  spectra,  depending  on  whether  longitudi¬ 
nal  or  lateral  correlations  appear  in  their  definition.  In  addition,  the  diagonal 
components  of  Eqs.  (20)  and  (21),  which  are  defined  by  autocorrelation  functions, 
are  usually  called  'power  spectral  densities'  or  'power  spectra'  while  the  off- 
diagonal  components  are  referred  to  as  'cross-spectral  densities'  or  simply  'cross- 
spectra'.  Since  the  spectra  of  the  0^.(n)  tensor  are  those  that  are  measured  in 

practice,  they  are  of  the  most  interest  and  unless  specified  otherwise,  the  terms 
'power “Spectra'  and  'cross-spectra'  refer  to  these  spectra  throughout  the  remainder 
of  this  review. 


Because  of  the  assumption  of  stationarity,  the  power  spectra  0i;i(n)  are 
even  functions  of  the  frequency  n.  Thus  from  Eqs.  (21)  and  (23), 


and 


/■>  00 

ft  (t)  cos(27nrr)dT 

o 

/■N  00 

0ii(n)  cos(27mT)  dn. 

n 


(25) 

(26) 


In  particular,  for  t  =  0,  one  obtains  from  Eq.  (26)  the  expressions 


&  (0)  s  u2  = 

'"■uu  ' 

r 00 

/  0  (n)  dn 

1  *uux  ' 

(a) 

“  o 

fe> 

o 

hi 

<rv>! 

it 

p  00 

I  tyyi*)  to 

(b) 

“  o 

&w<°)  ’  7  • 

f  *«,<">  I" 

(c) 

(27) 


indicating  that  the  area  under  the  power  spectrum  of  any  velocity  component  is 
equal  to  the  mean  square  of  the  component.  Equation  (25)  shows  that  the  power 
spectra  are  real  functions  of  n.  The  cross-spectra,  however,  are  complex  functions 
of  n,  with 

=  V(n)  sVn)* 1  + J- 


These  spectra  are  generally  redefined  as 


0ld(n)  =  £ 


ffhjO)  e 


■  i27mT 


dT 


=  C^(n)  -  i  ^j(n),  i  i  j 

where  C.^.(.n)  is  called  the  co-spectral  density  or  co-spectrum  and  Q^(n)  called 
the  qua  cl- spectrum,  and  these  functions  are  real-valued  even  and  odd1  functions  of 
n,  respectively.  In  addition,  these  cross-spectra  can  be  expressed  in  coefficient 
form  by  the  coherence  function,  defined  by 

?  I  ^(n)  I2 

hj(a)  "  OnC(n)  ’  1  t* 


where  y .  is  the  coherence. 
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Finally,  it  should  be  mentioned  that  in  the  spectral  theory  of  turbu¬ 
lence,  the  independent  frequency  variable  n  is  not  always  used.  In  the  litera¬ 
ture,  any  of  the  variables 

k  =  n/U,  w  =  27m,  or  fl  =  27m/U  =  277k 


may  be  found  where,  of  course,  if  Taylor's  Hypothesis  is  assumed,  k  3  k  .  In  these 
cases,  it  is  generally  Eq.  (2b)  that  is  used  to  properly  relate  the  various  spectral 

functions  to  each  other.  That  is 


r 


0ij(n)dn 


J  —  Q 


dk 


0^  (u>)  dw 


-*r 

J  -00 


dn 


which,  along  with  the  above  definitions  of  the  independent  variables,  can  be  used 

to  show  that 


dh  j ( k)  =  U  0  (n)  =  2ttU  ^(w)  =  2710  "(fl) 


(29) 


In  the  description  of  the  planetary  layer  which  follows,  it  is  the  function  0^ 
which  will  be  used.  This  is  because  if  the  'reduced  frequency'  or  inverse  wave¬ 
length  k  =  n/U  is  used  as  the  independent  variable,  where  0  is  the  mean  speed 
relative  to  the  measuring  system,  the  turbulence  measured  will  not  depend  on  the 
motion  of  this  system,  and  spectra  measured  by  aircraft  may  be  directly  compared 
with  those  obtained  on  towers  (at  the  same  altitude,  of  course).  Note  also  from 
Eq.  (29)  that 


00 


k  G^k)  =  n  0AJ(n) . 


2.5  Homogeneity  and  Isotropy 


The  turbulence  field  under  consideration  in  the  preceding  discussion 
was  assumed  to  be  homogeneous.  That  is,  its  statistical  properties  do  not  vary 
frcm  point  to  point  in  the  field,  and  thus  all  the  functions  described  are  inde¬ 
pendent  of  the  location  of  the  point  P  in  the  field.  In  addition,  the  assump¬ 
tion  of  homogeneity  allowed  the  correlations  of  Sec.  2.1  to  be  written  as  functions 
of  the  separation  between  the  points  rather  than  of  the  points  themselves,  and 
thus  simplified  the  description  considerably. 


The  concept  of  isotropy  simplifies  the  description  of  the  turbulence 
even  further.  If  a  turbulence  field  is  isotropic,  its  statistical  properties  are 
independent  of  direction  in  the  field,  and  thus  they  do  not  change  with  a  rota¬ 
tion  of  the  coordinates  axes.  Thus  isotropy  implies  homogeneity,  but  not  vice- 
versa.  All  off-diagonal  components  of  the  special  correlation  tensors  of  Sec. 

2.1  are  zero,  as  are  the  corresponding  cross-spectra  of  Sec.  2.b  (note  that  the 
off-diagonal  components  of  the  general  correlation  tensor  of  Eq.(^)  are  not 
necessarily  zero).  All  longitudinal  correlations  are  equal,  such  that 

f  (r  )  =  f  (r  )  =  f  (r  )  3  f(r)  (30) 

u  x  v  y  w  z7 

and  ail  lateral  correlations  are  equal,  or 

^v(rx}  r  gw(rx}  "gu(ry)  =*w'ry>  =  f5u(rz)  =  gv(rz)  3  g(r)  (3l) 
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and  the  equation  of  continuity  can  be  used  to  relate  these  correlations  by 

g(r)  *=  f(r)  +  |  (32) 

where  r  *  {  r|.  Thus  only  one  independent  non-zero  correlation  remains.  In 
addition,  Batchelor  (Ref. 7)  has  shown  that  in  isotropic  turbulence,  the  general 
(spatial)  correlation  tensor  is  given  in  terms  of  f(r)  and  g(r)  by 

r  r 

\j(£)  =  [f(*)-g(r)3  +  g(0  ,  i,J  =  u,  v,  w 

where  is  the  Kronecker  delta  and  ru  £  rx»  etc.  Similarly,  the  three- 
dimensional  spectrum  tensor  E^k)  nay  be  expressed  as  a  fairly  compact  function 
in  isotropic  turbulence  by 


where 


E  (*>. 
le-r'X 


and  E(A)  is  a  scalar  function  usually  referred  to  as  the  'energy  spectrum  function', 
for  the  0.,  (n)  spectra,  isotropy  requires  that 

0itJ(n)  =  0  ,  i  /  j 

and  thus  from  Eq.  (28), 

7i(j  *  0,  i  f  j  .  (33) 

Also, 

*w<n>  *  *w<n> 


and  Eq.(32)  may  be  used  to  show  that 

d 0  (n) 

0ww(n)  -  k  *, ,»  -  §  UU 


uu' 


dn 


(3*0 


Because  of  Eq:'.  (30}  and  (31),  all  longitudinal  integral  scales  are  equal  to  each 
other  in  isotropic  turbulence,  as  are  lateral  scales.  That  is, 


(35) 


ard 


ard  from  Eq.  (32^ , 


it  can  be  shown  that 


L  =  2  L  . 
u  v 


(36) 
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Finally,  the  velocity  component  mean-square  values  are  equal  in  isotropic  turbulence, 
such  that 


2 

u 


cj-2 


and  the  Reynolds  Stress  tensor,  for  example,  is  now  given  by 


(37) 


III,  CHARACTERISTICS  OF  THE  PLANETARY  BOUNDARY  LAYER 


3.1  General  Description 

The  atmosphere  near  the  surface  of  the  earth  can  conveniently  be  divided 
into  three  regions.  These  regions  are  the  free  atmosphere,  the  planetary  boundary 
layer,  and  the  surface  boundary  layer.  In  the  free  atmosphere,  viscosity  is  neg¬ 
lected,  and  only  inertial,  Coriolis,  and  pressure  gradient  forces  act  on  the  air. 

The  wind  resulting  from  these  forces  is  called  the  gradient  wind,  and  is  independent 
of  the  nature  of  the  earth's  surface  below.  It  can  easily  be  shown  (Ref. 3)  that 
the  gradient  wind  must  flow  along  the  isobars,  and  in  the  special  case  when  the 
isobars  are  straight  or  so  slightly  curved  that  centripetal  acceleration  of  the  air 
is  negligible,  the  gradient  wind  is  called  the  'geostrophic '  wind  and  is  given 
approximately  by 


where  n  is  normal  to  the  isobars  and  f  =  2wsin0,  with  w  the  rotational  velocity  of 
the  earth  and  0  the  latitude. 

The  planetary  boundary  layer  refers  to  the  region  between  the  earth's 
surface  and  the  height  at  which  the  free  atmosphere  can  be  said  to  begin.  This 
height  is  called  the  'gradient  height',  zQ  and  is  generally  of  the  order  of  1000- 

2000  feet,  depending  on  surface  conditions.  The  surface  boundary  layer  is  a  sub¬ 
layer  of  tne  planetary  layer,  extending  from  the  ground  up  to  about  200  ft.  (  + 

100  ft,)  depending  again  on  surface  conditions.  In  the  surface  layer,  Coriolis 
forces  are  assumed  negligible,  and  wind  characteristics  are  determined  by  surface 
conditions,  thermal  stability  and  height.  Shear  stress  is  assumed  constant  here, 
and  indeed  this  actually  defines  the  extent  of  the  region.  In  the  remainder  of 
the  planetary  layer  above  the  surface  layer,  Coriolis  forces  begin  to  have  an 
effect  on  the  wind,  surface  roughness  effects  decrease,  and  the  shear  stress  de¬ 
creases  from  its  constant  value  in  the  surface  layer. 

In  this  review,  it  is  the  planetary  layer  that  is  being  considered,  and 
consequently  the  basic  parameters  of  interest  will  be  the  surface  conditions, 
thermal  stability,  and  the  height. 


12 


) 


■  • Atmospheric  Stability 

The  hydrostatic  stability  of  the  atmosphere  depends  of  course  on  the  temp¬ 
erature  gradients  in  it.  The  basic  criterion  is  whether  or  not  the  decrease  in 
temperature  with  height,  called  the  lapse  rate  7 ,  is  greater  than  the  dry  adiabatic 
.apse  rate,  r(7  l°C/lOC  m.  5.5°F/lOOO  ft).  During  a  'lapse1  period,  when 

y  5  -dT/d2  >  T 

a  volume  of  air  displaced  upward  will  experience  a  buoyant  force  upward  and  thus 
will  continue  to  ascend.  Thus  the  atmosphere  in  lapse  periods  is  classified  as 
'hydrostatically  unstable',  or  simply  unstable.  During  so-called  'inversion' 
periods,  when  7  <  P,  the  volume  of  air  displaced  upward  will  be  at  a  lower  tern* 
peraiure  than  its  surroundings  and  will  experience  a  restoring  downward  force,  in 
which  case  the  atmosphere  is  classified  stable.  If  7  *  r  or  very  slightly  less, 
the  atmosphere  is  classified  neutral.. 

As  for  the  effect  of  stability  on  turbulence  in  the  atmosphere,  one  would 
obviously  expect  greater  turbulence  levels  in  hydrostatically  unstable  air  since 
under  these  conditions,  heat  convection  would  be  added  to  the  'mechanical'  turbu¬ 
lence  produced  by  the  shear.  A  better  understanding  is  obtained  if  one  considers 
the  driving  force  in  moving  vertically  the  volumes  of  air  mentioned  above.  The 
energy  required  to  displace  the  air,  that  is,  to  produce  turbulence,  is  extracted 
from  the  mean  flow  by  the  Reynolds  Stress.  Whether  the  velocity  fluctuations 
increase  or  decrease  will  depend  on  whether  or  not  the  rate  of  this  energy  supply 
is  greater  than  the  rate  at  which  work  must  be  done  in  the  gravitational  field 
r.u  moving  the  fluid  volumes  in  the  vertical  direction.  The  parameter  expressing 
this  criterion  is  the  gradient  Richardson  Number,  Ri,  defined  by 

Ri  *  .  (38) 

e(du/dz) 

Between  zero  and  unity  there  is  a  'critical'  Richardson  Number  above  which  turbu¬ 
lent.  motion  in  the  air  will,  subside  into  laminar  motion  and  below  which  it  will 
remain  turbulent.  4  definite  value  of  this  critical  number  is  not  available, 
and  indeed  it  may  depend  on  surface  conditions.  There  is  some  indication  (Ref's. 

3  and  4),  however,  that  it  should  be  Ri  ~  0.25,  such  that  above  this  value  of 
Ri  turbulence  ceases  to  exist. 

Hydrostatic  stability  can  easily  be  expressed  in  terms  of  Richardson 
Number.  For  an  ideal  gas,  the  potential  temperature  0  is  related  to  the  actual 
tempe  nature  T  by 


Thus  from  Eq.(38) 


1  be 

0  3z 


dT 

3? 


+  r 


Ri  - 


g(dT/dz  +  T) 
T(du/dz)? 

&iLz_2) _ 

T(dU/dz)? 


here fore  unstable  air  corresponds  to  <  0  and  stable  air  to  R  >  0. 


(39) 
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Summarizing,  then,  the  atmosphere  may  be  characterized  by  the  following  Richardson 
number  regimes:: 

Ri  <  0  unstable  air,  with  considerable  convective  turbulence  in 

addition  to  the  mechanical  turbulence 

R .  ~  0  generally  for  |Ri  0.03?,  the  air  is  termed  neutral  or 

' near- neutral'  and  the  turbulence  is  purely  mechanical 

(X  Ri  <~  0.25  stable  air,  with  mechanical  turbulence  being  damped  by 
the  thermal  stratification 

R.  >  ~  0.25  very  stable  air  in  which  no  turbulence  can  exist  at  all, 
at  least  in  the  vertical  direction 

The  Richardson  number  is  v:ry  important  in  the  atmosphere  since  the  Reynolds  numbers 
are  so  large  that  they  cease  to  be  of  importance,  and  Ri  is  the  most  relevant  non- 
dimensional  number. 

3.3  Simplifying  Assumptions 
3.3-1  Stationarity 

Spectral  data  measured  by  aircraft  at  different  times  in  the  planetary 
layer  nave  shown  that  stationarity  of  the  velocity  signals  can  reasonably  be 
assumed  for  periods  up  to  10-20  minutes,  and  sometimes  longer.  Thus  reasonable 
record  lengths  can  be  obtained  for  the  data  without  the  stationarity  assumption 
breaking  down. 

3.3*2  Homogeneity 

It  is  generally  assumed  that  the  flow  in  the  planetary  layer  is  homogeneous 
in  all  horizontal  planes,  such  that  aircraft  measurements  obtained  in  any  direction 
in  these  planes  may  be  treated.  Gunter  et  al  (Ref .43)  go  so  far  as  to  conclude 
complete  homogeneity  of  the  turbulence  based  on  measurements  at  250  ft. and  750  ft. 
However,  this  conclusion  seems  somewhat  optimistic,  and  indeed,  some  of  the  data 
presented  refutes  this  conclusion  since  they  find  a  definite  increase  in  integral 
scale  with  height.  Consequently  only  the  concept  of  horizontal  homogeneity  is  considered 
acceptable,  and  results  show  that  in  general  this  is  a  quite  reasonable  assumption, 
particularly  over  relatively  homogeneous  terrain. 

3.3.3  Isotropy 

In  the  free  atmosphere,  isotropy  is  a  fairly  reasonable  assumption,  but 
this  proves  less  tenable  with  decreasing  altitude.  Lappe,  Davidson  and  Notess 
(Ref. 19)  concluded  from  tower  and  aircraft  measurements  at  z  ~  300  ft. in  unstable 
air  that  Taylor's  Hypothesis  is  roughly  equally  well  satisfied  regardless  of  what 
direction  the  aircraft  flies  relative  to  the  wind.  This  would  indicate  that  the 
turbulence  under  these  conditions  is  more  or  less  horizontally  isotropic.  Gunter 
et  al  (Ref. 43)  tested  isotropy  for  hundreds  of  hours  of  aircraft  data  taken  at 
z  -  250  ft. and  750  ft.(i.e.,  generally  above  the  surface  layer)  by  comparing 
experimental  ratios  of  vertical-to- lateral  and  longitudinal-to- lateral  component 
spectra  with  the  corresponding  ratios  obtained  from  the  isotropic  von  K&rm6n 
spectral  equations.  Their  conclusion  is  that  for  most  stability  cases  and  for  all 
combinations  of  height  and  surface  conditions,  the  turbulence  is  totally  isotropic. 


2  2 

This  conclusion  is  enhanced  ty  y  y  ar.d  y  coherence  measurements,  which  were 

uv  vw 

always  less  than  0.1?  for  all  frequencies,  and  would  of  course  be  zero  for  iso¬ 
tropic  turbulence.  It  is  to  be  noted,  however,  that  the  7^  coherences,  which 
one  would  expect  to  be  larger  than  the  other  two  in  a  shear  flow,  were  not  pre¬ 
sented.  Also,  upon  inspecting  the  presented  data,  it  is  seen  that  the  conclusion 
of  isotropy  depends  largely  on  the  size  of  the  difference  between  the  theoretically 
isotropic  case  and  the  experimental  data  that  is  allowed  for  the  data  to  be  called 
isotropic.  Within  roughly  20$  variation,  the  data  is  indeed  seen  to  indicate  iso¬ 
tropy,  although  the  variation  of  scale  with  height  still  precludes  true  isotropy. 

It  it.  evident  from  the  data,  however,  that  there  is  a  distinct  reduced  frequency 
k  «  3  x  10" 3  cpf.  above  which  tuere  is  a  very  significant  decrease  in  the  amount 
of  departure  of  the  data  from  the  isotropic  case.  This  suggests  a  'local  isotropy' 
region  in  which  the  turbulence  is  significantly  closer  to  true  isotropy  than  at 
lower  frequencies. 

The  concept  of  local  isotropy  was  first  put  forth  by  Kolmogoroff.  It 
postulates  that  in  the  so-called  inertial  subrange  region  of  the  energy  spectrum, 
the  turbulence  is  isotropic.  In  this  region  turbulent  energy  is  neither  produced 
nor  dissipated,  but  merely  passed  through  from  the  large  anisotropic  eddies  to 
smaller  eddies  by  inertial  forces.  This  energy  is  then  dissipated  by  viscous 
forces  at  the  same  rate  e  at  which  it  is  inertially  passed  through  the  subrange, 
maintaining  an  equilibrium  state.  Kolmogoroff  showed  that  in  this  inertial 
subrange,  the  energy  spectrum  will  be  proportional  to  the  -?/3  power  of  frequency, 
and  the  _  tio  of  longitudinal-to- lateral  or  longitudinal- to- vertical  spectra  must 
be  3/4.  The  longitudinal  component  power  spectrum  in  this  region  is  of  the  form 

<t  (k)  -  b  e2/3  k“5/3 

where  b  is  a  constant,  while  the  lateral  and  vertical  component  spectra  are  given 
by  similar  expressions  but  with  b  replaced  by  4b/3.  In  the  planetary  layer,  this 
notion  of  local  isotropy  seems  quite  reasonable,  in  that  one  would  expect  anisotropy 
for  eddies  large  enough  to  be  affected  by  the  mean  shear  and  thermal  structure, 
while  smaller  eddies  with  shorter  time  scales  should  be  able  to  redistribute  their 
energy  among  the  components  more  quickly.  Thus  one  would  expect  that  over  a 
range  of  wavelengths  small  compared  with  some  characteristic  length,  say  the  dis¬ 
tance  to  the  ground,  z,  or  to  the  nearest  stable  layer,  the  turbulence  should  be 
isotropic . 


Experimental  results  have  in  general  tended  to  confirm  the  local  isotropy 
concept,  with  some  reservations.  The  spectra  generally  obey  the  -5/3  power  law, 
with  the  constant  b  having  a  value  of  ~  0.5  if  k  is  in  radians/m  or~0.06'5  if  k,is 
in  cycles/ ft.  However,  there  is  some  disagreement  as  to  whether  the  3/4  ratio  of 
the  spectra  in  the  subrange  actually  exists.  Elderkin  (Ref. 17)  found  for  spectra 
at  z  =  10  ft .  and  20  ft.  that  the  longitudinal  and  vertical  spectra  are  about 
the  same  in  magnitude  in  the  inertial  range,  and  refers  to  some  results  by 
Stewart  that  found  the  same  tendency.  Also  Berman  (Ref.40)  quotes  several  papers 
and  states  that  the  conclusion  from  these  results  is  that  all  component  spectra 
have  the  same  magnitude  m  the  subrange  and  not  the  Kolmogoroff  ratios .  On  the 
other  hand,  however,  Busch  and  Panofsky  (Ref. 41)  studied  considerable  data  for 
heights  up  to  300  ft.  and  concluded  that  in  regions  over  which  the  spectra  obey 
-5/ 3  power  laws,  the  ratios  of  the  components  "show  fair  agreement"  with  the  3/4 
ratios  predicted  by  Kolmogoroff  In  addition,  Fichtl  and  McVehil  (Ref. 39)  assumed 
the  3/4  ratio  m  obtaining  model  equations  for  the  longitudinal  and  lateral 
components  which  ultimately  fit  measured  data  very  well,  thus  indicating  this  ratio 
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to  be  valid  for  at  least  these  components.  Finally,  the  aircraft  data  obtained  by 
Gunter  et  al  obeyed  the  Kolmogoroff  ratios  quite  well  in  the  inertial  subrange.  In 
view  of  these  discrepancies,  no  definite  conclusion  can  be  drawn  as  to  whether  the 
Kolmogoroff  ratios  are  or  are  not  obeyed. 

There  is  far  better  agreement  among  investigators  as  to  the  validity  of 
the  -5/3  power  law  for  the  spectra.  The  data  generally  show  that  the  spectra  obey 
a  -5/3  frequency  dependence  up  to  wavelengths  much  greater  than  can  possibly  be 
expected  to  lie  in  the  Kolmogoroff  inertial  subrange,  especially  in  the  case  of  the 
longitudinal  spectrum.  However,  it  is  not  necessarily  the  case  that  this  represents 
a  simple  extrapolation  of  this  subrange  (see  Ref. 1*1).  As  for  the  actual  'isotropic 
limit'  above  which  local  isotropy  exists,  considerable  information  is  available. 
Lumley  and  Panofsky  (Ref.l*)  give  a  detailed  discussion  of  the  limiting  frequency, 
and  quote  a  result  from  Priestley  based  on  co-spectral  measurements  that  local 
isotropy  exists  for  k  >~  0.6/z.  Elderkin  (Ref .17)  concludes  that  while  the  -5/3 
law  is  obeyed  to  frequencies  as  low  as  k  ~0.2/z  in  some  cases,  the  co-spectra  do 
not  reduce  to  zero  until  k  ~3/z,  and  true  local  isotropy  is  not  attained  for 
k<  ~  U  /  z .  Lappe  and  Davidson  (Ref. 38)  state  that  the  Kolmogoroff  range  can  exist, 
if  at  all,  only  for  wavelengths  less  than  88  ft.  based  on  measurements  made  at 
z  =  1*00  ft.,  which  suggests  k  must  be  greater  than  ~ U/z .  Finally,  the  value  of  k 
quoted  above  from  the  da;ta  of  Gunter  et  al,  above  which  isotropy  improves  (k«3  x  10--3 
cpf.) ,  corresponds  to  an  isotropic  limit  of  k  ~l/z  -  2/z  for  the  heights  measured, 
a  value  roughly  in  agreement  with  that  suggested  by  Panofsky  in  Ref. 63.  Thus  it 
appears  that  a  reasonable  value  of  reduced  frequency  above  which  local  isotropy 
can  be  said  to  exist  is  k  ~3/z. 


To  summarize,  it  can  generally  be  said  that  turbulence  in  the  planetary 
layer  is  not  isotropic  in  the  true  sense.  However,  a  region  of  local  isotropy  does 
exist  for  k>~3/z,  although  the  component  spectra  may  or  may  not  obey  the  Kolmogoroff 
ratios  in  this  region.  Horizontal  isotropy,  in  which  the  turbulence  is  independent 
of  rotations  of  the  coordinate  system  about  the  vertical  axis,  is  also  a  reasonable 
assumption . 


3.3.1*  Taylor's  Hypothesis 

Lin  (Ref. 55)  has  investigated  the  validity  of  applying  Taylor's  Hypo¬ 
thesis  to  the  turbulence  in  shear  flows  and  concludes  that  the  requirement  of 
u'<<  0  should  be  valid  for  k  such  that 


dU 
dz  * 


(1*0) 


If  the  logarithmic  velocity  law  is  assumed  to  hold  (Eq.(l*2),  Sec. 3.1*)  then  Eq. 
(1*0)  becomes 


k  » 


1 

z  In  ( z/z  ) 
o 


(1*1) 


Lappe  and  Davidson  (Ref. 38)  compared  aircraft  and  tower  measurements  at  z  =  300-1*00 
ft.  and  found  that  the  spectra  so  measured  were  the  same  for  wavelengths  at  least 
as  large  as  600-900  ft.  Since  z  for  these  tests  was  ~3.3  ft.  (Ref. 19),  Eq.  (1+1) 
suggests  that  in  this  case,  k»~0.0005  rpf.  for  Taylor's  Hypothesis  to  be  valid, 
or 


A 


1 

k 


«  2000  ft. 
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Thus  Lappe 
it  must  be 
A  >  ~  1000 


and  Davidson's  result  roughly  obeys  Lin's 
realized  that  Taylor's  Hypothesis  becomes 
ft. 


requirement.  In  any  case, 
decreasingly  accurate  for 


3.4  Mean  Wind  Character is ties 


3.4.1  Velocity  Profiles 


In  the  surface  boundary  layer,  the  validity  of  the  Prandtl  logarithmic 
law  for  the  mean  velocity  in  a  neutral  atmosphere  has  been  well-verified.  The 
law  can  be  easily  obtained  from  Prandtl* s  Mixing  Length  theory  (or  von  K&rmfin's 
similarity  hypothesis)  under  the  assumptions  that 

(i)  viscous  stress  is  negligible; 


(ii)  mixing  length  is  proportional  to  height;  that  is,  l  =  Kz; 


(iii)shearing  stress  is  constant  and  equal  to  the  surface  stress  t  .  That  is, 


t  =  -puw  =  t  =  IT  p  _ 

O  T  r 

where  0  is  the  friction  velocity  and  is  defined  by  the  above  relation. 

If  in  addition  U  =  0  at  z  =  z  and  it  is  assumed  that  z  «  z,  the  law  can  be 
written  in  the  form 


where  K  is  the  von  Karmen  constant  (~0.4)  and  z  is  the  so-called  'roughness 
length' . 


The  logarithmic  law  of  Eq.(42)  is  valid  strictly  speaking  only  for 
neutral  conditions.  In  non-neutral  stabilities,  a  modification  of  this  law  is 
given  (Ref. 4  0"  44)  such  that 


U 


U 

T 


1 

K 


(43) 


where  L'  is  a  temperature  dependent  scaling  length.  The  function  ^  is  a  universal 
function  of  z/L*  and  since  it  can  be  shown  that  the  Richardson  Number  is  a  unique 
function  of  z/L'  (Ref .4),  Eq.(43)  can  be  used  to  obtain  mean  velocity  in  non¬ 
neutral  stabilities. 


The  roughness  length  z  usually  turns  out  to  be  about  l/30  times  the 
average  dimension  of  a  typical  rBughness  particle.  In  practice,  both  zq  and 

are  determined  f^om  measurements  of  at  least  two  values  of  U  and  z  in  the  surface 
layer  and  solution  of  Eq,  (43) .  An  excellent  description  of  the  details  of  the 
procedure  used  to  determine  these  quantities  with  the  greatest  possible  accuracy 
is  given  in  Ref.  44. 

Above  the  surface  layer,  Coriolis  forces  increase,  surface  r  .ghness 
effects  decrease,  and  the  logarithmic  law  begins  to  depart  from  the  empirical 
data.  However,  treating  the  planetary  layer  as  a  whole,  a  power  law  velocity 
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profile  of  the  form 


u/u1  =  (z/Z]?  (44) 

is  found,  to  fit  the  uata  quite  well  in  neutral  conditions,  where  U..  is  the  velo¬ 
city  at  some  reference  height  z,  .  Davenport  (Ref's.  9  and.  10)  uses  the  gradient 
height  as  the  reference  height ,  such  that 

U/UQ  -  (z/zG)a.  (45) 

The  parameters  zG  and  a  depend  on  surface  conditions,  and  Davenport  (Ref. 42)  has 
surveyed  published  wind  profiles  at  sites  having  a  wide  range  of  zq  values  in 
order  to  determine  the  dependence  of  these  parameters  on  z  .  The  results  of 
this  survey  are  given  in  Fig.  1,  and  the  shape  of  the  profile  given  by  Eq.  (45) 
for  three  different  terrain  types  ts  ■  shown  in  Fig. 2.  The  gradient  height  z  must 
in  general  be  obtained  by  measurement  of  U  at  several  heights  and  an  estimate 
of  the  gradient  velocity  0^.  Thus  there  is  some  uncertainty  in  its  exact  value 
since  Ur  is  usually  obtained  from  isobaric  charts  whose  accuracy  depends  on  the 
spacingJbetween  the  meteorological  stations  from  which  the  pressure  readings  are 
obtained.  However  recent  estimates  of  both  zG  and  a  have  shown  good  agreement 
with  the  values  suggested  by  Fig.  1,  including  the  results  for  various  cities 
around  the  world  given  by  Davenport  in  Ref.  10.  It  must  be  remembered,  of  course, 
that  the  power  law  can  be  of  only  limited  usefulness  in  cities,  since  it  cannot 
be  expected  to  apply  to  the  mean  wind  speed  for  heights  considerably  below  typical 
obstructions .  In  cities  these  obstructions  may  be  up  to  700-800  ft.  high,  and 
the  wind  speed  is  lery  much  a  f ’unction  of  the  detailed  nature  of  the  structures. 
Harris  (Ref. 58)  states  that  it  is  not  certain  at  the  present  whether,  for  the 
treatment  of  winds  above  a  city,  it  is  more  correct  to  use  a  power  law  with  large 
0c  as  Davenport  suggests, or  to  use  a  lower  index  combined  with  an  upwards  dis¬ 
placement  of  the  0=0  reference  plane  to  take  account  of  the  obstruction  height. 
However,  it  would  appear  from  the  data  that  Davenport's  approach  is  quite  reason¬ 
able  in  those  regions  where  any  power  law  can  be  expected  to  apply. 

In  applying  Eq.  (45)  to  estimate  gust  loads  on  buildings,  it  is  the 
extreme  values  of  the  mean  wind  that  are  of  importance.  Thus  the  value,  used 
for  U  is  in  general  based  on  surface  measurements  of  extreme  wind-speeds,  and 
Davenport  (Ref’s.  9?  10  and  42^  outlines  in  detail  the  procedure  for  estimating 

extreme  values  of  0_  from  local  meteorological  data. 

G 

As  for  the  effect  of  stability  on  the  shape  of  the  power-  law  profile, 
little  data  is  available.  In  general,  however,  the  exponent  a  tends  to  become 
smaller  with  decreasing  stability.  That  is,  in  unstable  air,  the  profile  is 
generally  'fuller'. 

3.4.2  Mean  Wind  Direction 


It  is  not  only  the  magnitude  of  the  mean  wind  in  the  planetary  layer 
which  is  important  but  also  the  direction,  particularly  in  the  consideration  of 
the  wind  loading  of  structures.  In  the  free  atmosphere,  the  mean  wind  is  the 
gradient  wind,  and  it  flews  along  the  isobars.  However,  in  the  planetary  layer, 
the  presence  of  shear  stresses  in  combination  with  the  Coriolis  forces  causes 
a  systematic  deflection  of  the  mean  wind.  This  deflection  is  away  from  the  iso¬ 
bars  in  the  direction  of  decreasing  pressure  gradient,  and  is  such  that  in  the 
Northern  Hemisphere  the  wind  direction  rotates  clockwise  with  increasing  height. 
Several  theories  have  been  postulated  as  to  the  exact  variation  with  height  of 
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the  angle  a  between  the  mean  wind  and  the  isobars  (i.e.,  the  geostrophic  direction), 
but  they  all  require  assumptions  about  the  exchange  of  momentum  due  to  turbulence 
at  various  heights.  This  exchange  is  represented  by  the  eddy  viscosity  where 


if  *  . 

*  ao/sz 

The  simplest  theory  is  due  to  Ekman  and  is  outlined  by  Sutton  in  Ref. 3.  Ekman . 
assumed  IC.  constant  with  height  and  as  a  result  predicted  that  a  should  decrease 
from  45°  at  the  surface  to  zero  at  the  gradient  height.  This  work  has  led  to  the 
change  of  wind  direction  with  height  being  referred  to  as  the  'Ekman  Spiral'. 

Ekman' s  theory  is  somewhat  of  an  oversimplification,  however,  and  Sutton  describes 
more  sophisticated  theories.  One  approach  by  Kohler  assumes  =  K^z111  where  m 

depends  mainly  on  stability,  and  predicts  that  at  the  surface, a  lies  between  10° 
and  30°,  decreasing  of  course  to  zero  at  zQ.  Sutton's  theory  predicts  values  of  a 
ranging  from  31°  in  open  country  to  45°  in  city  centres.  In  general,  although 
little  data  is  available  over  urban  areas,  other  measurements  have  indicated  values 
of  a  at  the  surface  of  the  order  of  20°,  agreeing  roughly  with  the  theory.  However, 
Harris  in  Ref .58  has  suggested  that  very  little  of  the  decrease  in  a  with  height  occurs 
in  the  first  600  ft.  or  so.  His  results  in  high  winds  (i.e.,  neutral  stability) 
over  flat  terrain  for  z  up  to  600  ft.  show  no  systematic  deviation  in  wind  direction, 
and  he  states  that  other  results  have  tended  to  agree.  Harris  thus  suggests  that 
except  for  very  tall  structures  or  those  with  special  features  making  them  extremely 
sensitive  to  wind  direction,  it  should  be  reasonable  to  ignore  the  change  of  wind 
direction  with  height  in  strong  winds  over  all  types  of  terrain. 

3.5  Reynolds  Stresses 

3.5.1  Total  Kinetic  Energy  of  Turbulence 

The  parameters  governing  the  magnitudes  of  the  variances  u',  v' ,  and  w'_ 
of  the  velocity  components,  and  thus  the  total  kinetic  energy  of  the  turbulence,  e, 
are  mean  velocity,  surface  roughness,  height,  and  stability.  Similarity  theory 
predicts  that  for  neutral  air,  the  total  kinetic  energy  is  roughly  proportional 
to  02,  indicating  that  this  energy  should  thus  be  independent  of  height.  Thus 

from  Eq.  (42)  we  get 

5  ecu?  -  -fSL_  (46) 

T  lr.2(z/z0) 

for  the  surface  layer,  and  the  data  (viz.  Ref. 4)  show  that  in  general  the  variation 
of  e  with  both  velocity  and  height  indicated  by  Eq.  (46)  is  roughly  correct.  The 
effect  of  surface  conditions,  however,  is  not  adequately  represented  by  the  rough¬ 
ness  length  z  .  This  result  is  not  surprising  since  z  is  a  measure  only  of  the 
small  scale  features  of  the  surface  and  does  not  account  for  large  irregularities 
such  as  hills  and  mountains.  As  for  stability,  a  decrease  tends  generally  to 
increase  the  total  energy,  due  to  the  additional  convective  energy  in  unstable 
air.  This  effect  is  minimized  over  rough  terrain,  where  most  of  the  turbulence  is 
of  mechanical  origin. 
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3.5.2  Vertical  Component  Variance 


Similarity  theory  predicts  that  in  neutral  stability, 

W  =  A  UT  (47) 

where  A  is  a  constant  independent  of  height,  windspeed,  and  surface  roughness. 

The  data  of  various  investigators  (.viz.  Ref's.  17,  29,  and  4l)  has  shown  this  to 
be  true  at  least  through  the  surface  layer,  indicating  that  the  effect  of  surface 
conditions  on  w'  is  adequately  represented  by  Ut,  and  thus  by  zq  (from  Eq.  (42)). 

That  this  is  so  for  the  w'  component  likely  results  from  the  fact  that  very  large 
wavelength,  low  frequency  components  of  w'  are  suppressed  due  to  the  presence  of 
the  ground,  and  thus  the  effect  of  large  scale  non- uniformities  is  reduced. 

As  for  the  value  of  the  constant  A  in  Eq.  (47),  there  is  considerable 
discrepancy  in  the  data.  Lumley  and  Panofsky  (Ref. 4)  point  out  that  values  rang¬ 
ing  from  0.7  to  1.3  have  been  obtained,  and  suggest  that  A  ~  1.0  is  the  best  com¬ 
promise.  However,  more  recent  data  has  leaned  toward  a  value  of  ~  1.3.  Elderkin 
(Ref .17)  obtained  a  value  of  1.33  for  z  --  10-20  ft;,  while  Busch  and  Panofsky 
(P’.ef.4l)  obtained  1.29  from- an  integration  of  their  model  spectrum,  and  state  that 
this  value  is  consistent  with  direct  measurements  made  using  sonic  anemometers. 
Also,  Panofsky  and  McCormack  (Ref .29)  found  a  value  of  ~  1.3  for  heights  up  to 
about  250  ft.  and  Panofsky  in  Ref .63  shows  some  results  indicating  that  A  =  1.3 
over  a  considerable  range  of  Richardson  numbers.  Thus  it  is  concluded  that 
A  =  1.3  is  the  most  appropriate  value  in  view  of  the  existing  data. 

The  effect  of  atmospheric  stability  on  w'  is  in  general  not  too  large. 
The  magnitude  of  w‘  tends  to  increase  somewhat  with  decreasing  stability,  as  one 
would  expect.  In  addition,  w’  increases  slowly  with  height  in  unstable  air, 
while  in  stable  air  it  decreases  with  height.  Thus  stability  effects  on  w'  tend 
to  be  felt  more  at  higher  values  of  z,  and  are  not  too  important  near  the  ground. 

3-5.3  Lateral  Component  Variance 

In  neutral  air,  at  least,  most  of  the  existing  data  show  that  through 

the  surface  layer, 


V  =  B  UT  (48) 

where  B  is  a  constant.  This  s;iggests  that  v'  does  not  vary  with  height  in  neutral 
air.  One  notable  exception  to  this  conclusion  is  the  investigation  by  Fichtl  and 
McVehil  (Ref. 39)  whose  tower  data  measured  for  heights  up  to  ~  500  ft.  indicate 
that  in  neutral  air,  v'  decreases  slightly  with  height,  being  proportional  to 
z-.17o  over  this  range.  As  for  the  dependence  of  v'  on  surface  conditions,  the 
roughness  z  is  unfortunately  inadequate  to  completely  estimate  the  surface 
effect,  as  is  the  case  for  the  total  kinetic  energy.  Consequently,  the  constant 
B  is  a  function  of  large  scale  roughness  and  thus  varies  from  place  to  place, 
ranging  from  1.3  to  2.6  (Ref  .4). 

Atmospheric  stability  has  a  quite  large  effect  on  v' ,  increasing  it  by 
a  large  amount  as  stability  decreases  for  the  same  wind  speed.  Lumley  and 
Panofsky  (Ref. 4)  show,  however,  that  in  either  stable  or  unstable  air,  there  is 
still  very  little  vertical  variation  of  v' .  In  this  case,  Fichtl  and  McVehil's 
results  agree  since  they  find  that  in  unstable  air,  v' oc  z-1^*®2,  indicating  a 
very  weak  dependence  of  v'  on  height. 
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An  interesting  consequence  of  the  differential  behaviour  of  v1  and  w' 

with  height  in  stable  air  is  pointed  out  by  Lumley  and  Panofsky.  Since  w'  decrease* 
with  height  in  stable  air  while  v'  does  not,  smoke  emitted  at  z  ~  300  ft.  in  an 
inversion  meanders  horizontally,  with  little  vertical  spreading.  In  fact,  this 
meandering  results  from  the  presence  of  low  frequency  components  in  the  spectrum 
of  v  (see  Sec.  3*6.2), 

3.5.4  Longitudinal  Component  Variance 

While  the  longitudinal  fluctuating  velocity  component  is  the  largest  of 
the  three,  its  properties  generally  tend  to  be  intermediate  between  those  of  v 
and  w.  Stability  affects  it  more  than  the  w  component  but  not  so  drastically  as 
it  affects  the  v  component.  Lumley  and  Panofsky  (Ref .4)  indicate  that  u'  does 
not  vary  with  height  ir.  any  stability,  although  the  fluctuations  in  stable  air 
tend  to  shift  from  high  frequency  near  the  ground  to  lower  frequency  at  larger 
heights.  The  independence  of  u'  with  height  is  confirmed  by  data  quoted  by  Templin 
(Ref .30)  for  unstable  air  and  heights  up  to  500  ft.  and  by  the  results  of  Fichtl 
and  McVehil  (Ref. 39)  who  for  similar  conditions  find  that  u'oC  z-0*^.  In  neutral 
stability, however ,  Fichtl  and  McVehil  find  that  u'oc  z-0.315,  rather  than  being 
constant  with  height. 

In  neutral  air,  Lumley  and  Panofsky  suggest  that 

u'  =  C  UT  .  (49) 

As  with  the  lateral  component,  the  roughness  length  z  is  again  not  adequate  to 
completely  represent  large  scale  variations  in  the  terrain,  and  while  C  is  constant 
with  height,  it  varies  from  place  to  place.  Values  quoted  range  from  2.1  to  2.9, 
with  the  value  of  2.5  suggested  as  a  reasonable  engineering  approximation. 

3.5.5  Summary  of  Variances 


The  preceding  discussion  of  the  three  component  variances  referred  for 
the  most  part  to  the  surface  boundary  layer.  In  this  layer,  if  neutral  stability 
is  assumed  and  the  suggested  values  of  the  constants  A,  B,  and  C  are  used,  the 
ratios  of  the  component  variances  are  given  by 

u'/v'/w'/UT  =  2. 5/2. 0/1. 3/1  .  (50) 

Also,  since  the  logarithmic  law  of  Eq.(.42)  is  valid  under  these  conditions,  the 
friction  velocity  is  given  by 

n  =  Q _ 

T  2.5  In  (z/zq) 

and  thus  from  Eqs.  (47),  (48),  and  (49)  the  component  intensities  are  given  by 


u 5  a  1 

U  In  (z/z  ) 
o 


V  =  0.80 

U  ln(z/zQ) 


and 


w1 

U 


0.52 

ln(z/z  ) 
'  0 


(51) 


These  curves  are  shown  in  Fig. 3  for  z  =3  cm,  a  value  characteristic  of  flat, 
open  country.  Also  shewn  are  some  actual  measurements  quoted  in  Ref. 30  for 
Sale,  Australia  which  has  similar  surface  conditions  (z  =  1  cm),  and  these 
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values  are  seen  to  agree  well  with  the  predicted  values. 

For  the  planetary  layer  above  the  surface  layer,  little  information  has 
been  presented  because  there  does  not  appear  to  be  a  great  deal  of  information 
available.  This  may  result  from  the  fact  that  most  towers  used  for  measurement  do 
not  reach  heights  far  above  the  surface  layer,  so  aircraft  data  is  the  only  major 
source.  Borne  useful  measurements  have  been  obtained  by  Harris  (Ref. 65)  for  heights 
up  to  about  600  ft.  in  neutral  stability,  and  these  are  shown  in  Fig. 3.  It  is 
seen  that  these  results  agree  quite  well  with  Eq.(5l)  for  u’/U  in  the  surface  layer 
and  with  the  Sale  data  above  it. 

Considerable  data  on  the  component  variances  at  z  =  250  ft.  and  750  ft. 
is  to  be  found  in  Ref .43.  Unfortunately,  since  velocity  values  are  not  available, 
direct  comparison  with  the  above  data  cannot  be  made.  However,  the  data  suggests 
that  both  these  heights  be  considered  above  the  surface  layer,  and  if  this  is  done 
some  useful  conclusions  can  be  drawn  about  u' ,  v' ,  and  w'  above  the  surface  layer. 
First,  within  reasonable  limits,  u'  «=  v,R=  w'  at  each  height.  Also,  increased  sur¬ 
face  roughness  results  in  larger  component  variances,  although  the  effect  is  felt 
more  at  250  ft.  than  at  750  ft., as  expected.  In  addition,  for  both  heights  there 
is  an  increase  in  the  variance  with  decreasing  stability,  and  finally,  the  variances 
decrease  from  250  ft.  to  750  ft.  for  virtually  all  stabilities.  When  these  con¬ 
clusions  are  considered  together  with  the  above  results,  it  would  appear  that  in 
general,  above  the  surface  layer,  the  three  velocity  component  variances  tend  to 
decrease  with  height  toward  the  same  value,  as  is  the  case  in  a  typical  wall 
boundary  layer  (viz.Ref.2  or  60) .  Thus  curves  similar  to  those  found  in  a  wall 
boundary  layer  have  been  shown  in  Fig. 3.  The  particular  value  that  the  variances 
approach  must  of  course  be  that  in  the  free  atmosphere,  which  is  very  small  and 
which  would  in  the  ideal  case  of  no  velocity  gradients  be  zero. 

3.5.6  Other  Reynolds  Stresses 

The  component  variances  discussed  above  are,  of  course,  singly  the  roots 
of  the  diagonal  components  of  the  Reynolds  Stress  tensor,  Eq.  (!?).  The  remaining 
Reynolds  Stresses  are  the  uv,  vw  and  uw  co-variances  (the  tensor  is  Hermitian), 
which  would  of  course  be  zero  if  the  turbulence  were  isotropic.  Elderkin  (Ref. 17) 
obtained  values  of  these  stresses  in  various  stabilities  at  z  =  10  ft.  and  20  ft. 
over  flat  terrain  having  zq  =  3  cm.  His  results  in  neutral  air  show  that  the  vw 

stresses  are  quite  small  compared  to  uw,  while  uv  values  are  larger  than  vw  but 
still  considerably  less  than  uw.  For  the  uw  stress,  there  is  really  not  sufficient 
data  to  make  significant  conclusions  about  the  assumption  of  constant  stress  in 
the  surface  layer,  although  the  values  at  the  two  measurement  heights  do  not  in 
general  show  any  larger  variation  than  that  between  different  runs  at  the  same 
height . 


In  the  surface  layer,  the  friction  velocity  is  given  by  (see  Sec.  3*4.1), 


U 


T 


Thus  if  the 
Eq.(50)  are 
is  given  by 


air  is  neutrally  stable  and  the  ratios  of  the  component  variances  of 
assumed,  then  the  correlation  coefficient,  for  the  u  and  w  components 
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u'  w' 


(2.5Ut)(1.3Ut) 


=  -  0.31. 


(52) 
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Above  the  surface  layer,  this  value  must  decrease,  since  the  shear  decreases  and 
uw  falls  to  zero.  As  for  uv  and  vw  above  the  surface  layer,  little  data  is  avail* 
able,  although  the  coherence  data  of  Gunter  et  al  (Sec,  3.3.3  and  3.6. 4$  suggest* 
that  these  stresses  are  quite  small. 

3.6  One-Dimensional  Spectral  Density  Functions 

In  general,  the  one-dimensional  turbulence  component  power  spectral 
densities  $^00  consist  of  three  distinct  regions:  a  low  frequency  region  in 
which  energy  is  put  into  the  turbulent  motion  from  the  mean  flow;  a  middle  fre¬ 
quency  equilibrium  region  in  which  no  energy  production  or  dissipation  is  occur¬ 
ring,  and  only  an  inertial  transfer  of  energy  toward  smaller  eddies  is  taking 
place  (i.e.,  the  Kolmogoroff  inertial  subrange);  and  a  high  frequency  region  in 
which  turbulent  energy  is  converted  to  heat  by  viscous  dissipation  (viscous 
subrange).  The  frequency  dependence  of  ^^(k)  in  these  three  regions  is  roughly 

given  by  k°,  k"5/3,  ^-7^  respectively.  As  with  the  component  variances, 

the  power  spectra  in  general  are  dependent  on  height,  surf ace. roughness,  mean 
velocity,  and  thermal  stability. 

In  addition  to  the  above  general  regions  of  the  power  spectra,  the 
existence  of  a  'buoyant  subrange'  has  been  suggested  for  atmospheric  turbulence 
spectra  in  stable  conditions.  In  this  region,  &  'buoyant  length'  depending  on 
the  stratification  is  the  factor  limiting  the  inertial  subrange,  rather  than 
the  distance  to  the  ground.  Lumley  ana  Panofsky  (Ref .4)  discuss  the  buoyant 
subrange  at  some  length,  and  Busch  and  Panofsky  (Ref.4l)  show  experimental  data 
suggesting  that  it  does  exist,  at  least  for  the  vertical  spectra,  with  a  frequency 
dependence  of  the  order  of  k"3  as  suggested  by  Lumley  (Ref. 56). 


The  behaviour  of  the  power  spectra  in  the  inertial  subrange  has  been 
discussed  at  some  length  in  Section  3.3.3  with  regard  to  both  the  relative  mag¬ 
nitudes  of  the  components  and  their  frequency  dependence,  and  consequently  there 
is  only  minimal  attention  given  to  this  region  in  the  discussion  that  follows. 


3.6.1  Vertical  Component  Spectrum 

The  vertical  velocity  power  spectrum  ^^(k)  has  been  quite  extensively 

investigated  and  a  considerable  amount  is  known  about  it.  In  general,  it  is 
found  that  below  a  height  of  roughly  150  ft. -200  ft.  the  spectrum  obeys  dynamic 
similarity.  That  is,  it  can  be  expressed  as 


k<t>  (k) 
wwv  _ 

=# 

W*" 


=  F(f,  Ri) 


(53) 


where  F  is  a  universal  function,  and  f  =  kz, 
plotted  in  similarity  coordinates  kO^k)/^ 


such  that  when  the  spectrum  is 
vs.  kz  for  fixed  stability,  it  is 


independent  of  height,  mean  velocity,  and  surface  roughness.  In  addition,  the 
shape  of  the  spectrum  at  these  heights  does  not  vary  a  great  deal  with  stability, 
except  perhaps  in  the  unstable  case  when  distinct  low  frequency  convective  peaks 
are  found  (viz* Ref. 17  for  z  =  40  ft.,  and  also  z  -  270  ft.).  The  maximum  value 
of  the  non-dimensional  spectrum  as  given  by  Eq.  (53)  was  found  to  be  ~  0.2  by 
Elderkin  and  slightly  larger  ( .23  -  .3)  by  Busch  and  Panofsky  (Ref.^l). 


While  stability  does  not  in  general  have  too  great  an  effect  on  the 
shape  of  the  spectrum  at  lew  altitudes,  it  does  tend  to  shift  the  spectrum  as  a 
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whole.  A  decrease  in  stability  moves  it  toward  lower  frequencies  (larger  wave¬ 
lengths)  while  an  increase • shifts  it  to  higher  frequencies.  This  is  a  reasonable 
effect,  in  that  in  stable  air,  convective  energy  is  minimized  and  one  would  thus 
expect  a  reduction  in  the  low  frequency  energy  components.  Thus  if  f  w=  k  w.  z 

is  the  value  of  f  at  which  the  non-dimensional  spectrum  has  its  maximum  vafue, 
the  effect  of  stability  is  indicated  by  changes  in  value  of  f  w.  Above  l60  ft., 
Busch  and  Panofsky  have  noted  that  stability  affects  not  only**  f  w,  but  also  the 
spectrum  shape,  with  energies  being  concentrated  around  higher  frequencies  in 
stable  air. 


For  regions  close  to  the  ground,  one  would  expect  the  energy  at  low 
wave  numbers  to  be  limited  by  the  presence  of  the  surface.  This  suggests  that 
the  scale  of  the  vertical  component  should  increase  with  height.  This  is  indeed 
the  case  since,  as  indicated  by  the  similarity  of  the  spectrum,  f  w  does  not 

change  with  height  for  fixed  stability.  Thus  k^w  is  inversely  proportional  to 

height,  and  since  it  is  also  inversely  proportional  to  scale  (see  Sec. 3. 7),  the 
scale  is  linearly  dependent  on  height.  Above  160  ft.,  Busch  and  Panofsky  have 
shown  that  f^w  begins  to  increase  with  height,  indicating  the  breakdown  of 

dynamic  similarity  and  a  weakening  dependence  of  scale  on  height.  More  will  be 
said  concerning  the  height  dependence  of  scale  in  Sec.  3*7* 


Many  analytical  expressions  have  been  suggested  in  the  past  in  order 
to  represent  the  power  spectrum  of  the  vertical  component.  Bowne  and  Anderson 
(Ref. 5)  list  several  of  these  expressions,  and  in  their  Fig.  15  compare  them 
for  similar  conditions.  Lumley  and  Panofsky  (Ref .4)  also  discuss  some  of  these 
spectra  in  detail.  One  well-known  model  is  the  Dryden  spectrum,  which  can  be 
obtained  for  isotropic  turbulence  by  assuming  that  the  correlation  function 
Rm(r  )  is  given  by 
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Then  using  Eq.(32),  the  correlation  R  (r  )  may  be  obtained  as 
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and  the  ^ww(kx^  Power  spectrum  is  found  by  Fourier  transforming  this  correlation. 
Assuming  Taylor's  Hypothesis,  this  spectrum  is  given  in  the  frequency  domain 

by 
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where  L  is  the  longitudinal  integral  scale  and  has  been  assumed  to  be  twice 
the  vertical  component  scale  L  x.  This  model,  like  several  others,  indicates 
a  frequency  dependence  of  k“?  ir  the  inertial  subrange,  and  because  the  recent 
data  consistently  show  a  -5/3  power  law,  these  models  have  fallen  out  of  use. 

Of  the  other  models,  the  von  K&'-man  equation  is  likely  the  best,  since  it  captures 
the  features  of  scale  length,  total  variance,  and  the  -5/3  slope.  Like  the 
Dryden  model,  it  was  developed  for  isotropic  turbulence,  and  L  x  =  2  Lwx  is 

assumed.  The  equation  is 


2b 


Another  good  model  Is  that  suggested  by  Eusch  and  Panofsky  (Ref.4l)  for  z  up  to 
160  ft*,  nondimensionally  given  by 
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assuming  w'  =  1.3  U  .  In  this  expression,  f  w  is  a  function  of  stability  only, 
and  Busch  and  Panofsky  relate  fpW  and  stability  in  their  paper. 


For  heights  above  the  surface  layer,  the  results  presented  by  Gunter 
et  al  (Ref .43)  leave  little  doubt  that  the  von  K&rmfin  model  spectrum  best  fits 
the  experimental  data.  For  lower  heights,  Fig, 4  shows  that  while  the  Busch- 
Panofsky  model  of  Eq.  (57)  provides  a  somewhat  better  fit  to  the  data  than  the 
von  K£rm&n  equation  at  low  frequencies,  both  models  tend  to  underestimate  the 
energy  in  this  region.  Thus  in  the  interest  of  simplicity  and  since  the  scale 
parameter  in  Eq.  (56)  can  be  used  to  indicate  stability  effects  as  is  f  in  Eq. 
(57),  it  is  suggested  that  the  von  K£rm£n  model  with  an  appropriate  valfie  of 
scale  (see  Sec.  3*7  and  4.6)  adequately  describes  the  vertical  component  spec¬ 
trum  throughout  the  planetary  layer. 


3.6.2  Lateral  Component  Spectrum 

The  power  spectrum  (k)  of  the  lateral  velocity  component  does  not 

in  general  obey  similarity  theory,  except  for  large  wave  numbers  in  some  cases. 
Also,  in  neutral  and  unstable  conditions,  there  is  little  dependence  of  the  shape 
of  the. spectrum  on  height,  although  the  location  of  the  peak  k  v  does  appa¬ 
rently  depend  on  height  as  suggested  by  the  results  of  Fichtl  and  McVehil  (Ref. 
39)>  who  found  that  for  z  =  60  -500  ft.  in  neutral  air,  k  voc  z  -0.42  an(j  in 
V  o  P 

unstable  air,  k  oCz”u'  °.  There  is,  however,  a  very  large  variation  in  both 

the  magnitude  and  shape  of  this  spectrum  with  stability.  In  addition  to  a  ten¬ 
dency  for  the  spectrum  to  shift  as  a  whole  toward  lower  frequencies  for  un¬ 
stable  conditions,  decreasing  stability  also  greatly  increases  the  low  frequency 
portions  of  the  spectrum  while  leaving  the  high  frequency  part  relatively  un¬ 
affected.  Thus  for  fixed  wind  speed  the  effect  of  decreasing  the  stability 
is  to  superimpose  long  period  variations  on  the  shorter,  mechanically  produced 
eddies.  Consequently,  the  low  frequency  part  of  the  spectrum  depends  mainly 
on  stability  while  the  high  frequency  part  depends  on  roughness  and  velocity. 


In  stable  air,  height  variation  of  ^^(k)  becomes  more  significant. 

Lumley  and  Panofsky  (Ref.-:-)  show  data  which  indicates  that  while  the  low  frequency 
part  of  the  spectrum  remains  roughly  constant  with  height,  the  high  frequency 
part  decreases  rapidly.  The  relative  invariance  of  the  low  frequency  part  sug¬ 
gests  also  that  the  total  variance  v^  tends  not  to  vary  too  much  with  height,  as 
indicated  in  Sec.  3.5.3.  Recalling  that  w^  decreases  with  height  in  stable  air, 
and  since  k<t>ww(k)/w?  is  invariant  with  height  (dynamic  similarity)  it  follows 

that  $ww(k)  decreases  with  height  for  all  frequencies.  Thus  at  heights  of  the 
order  of  200-300  ft.  in  stable  air,  the  low  frequency  components  of  <t^(k)  are 


large  relative  to  those  of  <t  (k),  explaining  why  smoke  plumes  in  these  conditions 
meander  later ially  with  littYe  vertical  spreading. 


The  lateral  component  spectrum  is  more  difficult  to  represent  analytically 
than  that  of  the  vertical  component  because  of  the  lack  of  similarity.  The 
usual  practice  ±6  to  use  the  same  model  as  for  the  vertical  spectrum,  btt  with 
w2  replaced  by  v^  and  Lwx  replaced  by  V-  In  the  case  of  the  von  Kflrm&n  model, 

this  theoretically  makes  no  difference,  since  in  isotropic  turbulence, 
LX-2LX=2LX  and  =  w-, and  the  equation  is 

U  V  W 


$  (k)  =  2  v 

w'  ' 


-o  x  |1  +  188.4  (Lux  k)2 
u  H1+  70.7(Lux  k)2]11/6 


} 


(58) 


Another  equation  that  has  been  found  appropriate  is  a  modification  of  the  Busch- 
Panofsky  relation  of  Eq.(57)  which  Fichtl  and  McVehil  fitted  to  their  data. 

This  equation  is  given  non-dimensionally  by 


k  <t>  (k)  C  f/f  V 

w  _  v  '  p 

[1  +  i.5(f/f  v) 

where  and  ry  depend  on  stability  and  f^v  and  depend  on  height  and  stability. 

Above  the  surface  layer,  the  results  of  Gunter  et  al  (Ref .43)  show  that 
the  von  Kerman  equation  provides  the  best  fit  to  the  data  for  a  wide  range  of 
conditions.  Consequently  it  is  suggested  that,  as  with  the  vertical  component, 
the  von  Karirfan  equation  be  used  to  model  the  lateral  component  spectrum  through¬ 
out  the  planetary  layer,  using  appropriate  values  of  scale  (see  Sec. 3. 7  and  4.6). 

3.6.3  Longitudinal  Component  Spectrum 

3. 6. 3.1  Total  Velocity  Spectrum 

The  complete  power  spectrum  of  the  total  velocity  in  the  longitudinal 
direction  extends  to  much  lower  frequencies  than  those  of  the  other  two  velocity 
components.  This  is  due  to  the  long  period  variations  of  the  ’mean’  wind  speed 
0  resulting  from  macrometeorological  fluctuations,  and  U  is,  of  couae,  non¬ 
zero  only  in  the  longitudinal  direction.  The  best  data  available  for  the  complete 
spectrum  of  horizontal  windspeed  is  that  obtained  by  van  der  Hover*  (Ref. 32)  at 
Brookhaven,  and  shown  in  Fig,  5.  This  spectrum  was  formed  by  piecing  together 
several  separate  portions  which  were  obtained  at  different  times  at  z  ~  330  ft. 
Also  shown  in  Fig.  5  for  comparison  are  the  periods  of  some  common  occurrences 
as  obtained  from  Ref.  30, in  addition  to  a  pair  of  conjectural  peaks  suggested  by 
Davenport  (Ref. 9)  and  not  measured  by  van  der  Hoven. 

Generally,  the  part  of  the  complete  spectrum  measured  by  van  der  Hoven 
displays  two  major  peaks.  One  of  these  peaks  has  a  period  of  ~  100  hours,  and 
is  characteristic  of  the  large  scale  movement  of  pressure  systems.  The  other 
major  peaK  has  a  period  of  1-2  minutes  and  is  characteristic  of  the  gusts  and 
winds  resulting  from  the  modifying  influence  of  the  ground.  It  is  this  part  of 
the  curve,  the  ’gust1  spectrum,  which  is  represented  by  the  fluctuating  component 
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u(t)  and  which  corresponds  to  the  lateral  and  vertical  component  spectra  already 
discussed.  It  is  of  major  interest  since  it  contains  the  frequency  components 
which  most  seriously  affect  structures  and  aircraft.  In  between  the  two  peaks, 
there  is  a  minor  peak  with  a  period  of  about  twelve  hours  which  corresponds  to 
the  lulls  in  windspeed  generally  occurring  at  sunrise  and  sunset. 

Between  tne  gust  peak  and  the  twelve  hour  peak,  there  is  an  obvious 
region  of  very  row  energy  extending  from  a  period  of  ~5  minutes  to  ~5  hours. 

This  region  is  referred  to  as  the  'spectral'  or  'micrometeorological  gap'  and 
it  apparently  exists  for  all  stability  and  terrain  conditions.  It  results  from 
a  lack  of  any  physical  processes  capable  of  generating  fluctuations  having. these 
periods,  and  is  generally  assumed  to  exist  for  all  cases.  This  assumption  does 
not  take  into  account  the  occurrence  of  intense  local  storms  such  as  thunderstorms 
and  tornadoes,  but  these  are  usually  of  such  local  extent  that  their  probability 
of  occurrence  at  a  particular  location  is  negligibly  small. 

The  existence  of  the  spectral  gap  is  indeed  fortunate  in  that  it  allows 
a  very  clear  distinction  to  be  made  between  gusts  and  macrometeorological 
fluctuations.  By  defining  the  mean  wind  velocity  0  as  the  average  over  some 
period  within  the  gap,  it  is  conveniently  separated  from  the  gust  velocities. 

Also,  since  the  gap  begins  at  periods  as  short  as  5  minutes,  estimates  of  0  can 
be  obtained  from  very  reasonable  record  lengths.  Indeed  the  meteorological  esti¬ 
mates  of  mean  wind  in  various  countries  throughout  the  world  are  based  on  avera¬ 
ges  over  periods  ranging  from  5  minutes  to  one  hour,  and  the  existence  of  the 
spectral  gap  ensures  that  these  estimates  will  be  virtually  the  same. 

3.6.3.?  Gust  Velocity  Bpectrum 

The  spectrum  <t>  (k)  of  the  longitudinal  gust  component  in  the  surface 

layer  is  in  a  sense  intermediate  between  that  of  the  vertical  and  lateral  velo¬ 
city  components.  The  Tow  frequency  portion  is  affected  more  by  stability  changes 
than  the  high  frequency  region,  so  the  shape  of  the  spectrum  is  stability  depen¬ 
dent.  However,  the  dependence  is  not  nearly  as  strong  as  that  of  the  lateral 
component,  while  the  vertical  component  shape  is  virtually  independent  of  stabi¬ 
lity.  It  is  noted  also  that,  as  for  the  other  velocity  components,  decreasing 
stability  shifts  the  spectrum  to  lower  frequencies. 

Under  ail  conditions,  the  longitudinal  component  spectrum  contains 
considerable  energy  at  low  frequencies,  such  that  in  stable  and  neutral  air, 
when  this  energy  is  less  than  in  unstable  air,  it  is  still  considerably  larger 
than  that  of  the  other  two  component  spectra.  This  is  an  indication  that  in 
neutral  and  stable  air,  the  larger  eddies  are  elongated  in  the  wind  direction 
(see  also  Sec.  3*7).  Another  distinct  feature  of  the  longitudinal  component 
spectrum  is  that  the  -5/3  power  law  extends  to  considerably  larger  wavelengths 
than  it  does  for  the  other  two  spectra.  Lumley  and  Panofsky  (Ref.1*)  suggest 
it  is  obeyed  down  to  k  ~  0.2/z,  the  value  also  found  by  Elderkin  (Ref. 17)  and 
by  Panofsky  and  van  der  Hover.  (Ref. 35)* 

In  general,  as  with  $  (k),  the  spectrum  does  not  obey  dynamic  simi¬ 

larity  except,  perhaps,  at  very  large  wave  numbers.  When  plotted  in  similarity 
coordinates,  the  spectra  show  a  great  variety  of  sh^pe  and  position  for  neutral 
air  (viz.  Ref.4o).  The  maximum  value  of  k  <t>uu(k)/uf-  also  varies,  since  Busch 

and  Panofsky  (Ref ,4j )  found  a  value  of  roughly  0.18,  assuming  u'  -  2.5  0  ,  and 
Berman  (Ref.4f)  obtained  ~  0.27  using  the  same  a-  umption,  while  Elderkin  (Ref. 17) 


found  a  value  of  0.12  and  Fichtl  and  McVehil  (Ref .39) found  ~  0.13.  Finally,  the 
location  of  the  spectral  peak  f  u  violates  dynamic  similarity,  in  that  its  value 
tends  in  general  to  increase  wi?h  height  at  all  locations.  This  result  is  found 

by  virtually  all  investigators,  and  as  the  value  of  f  u  determines  the  scale  of 

P 

the  longitudinal  component,  its  variation  will  be  discussed  further  in  Sec.  3.7. 


As  has  been  done  for  the  vertical  component  spectrum,  a  number  of 
analytical  expressions  have  been  suggested  to  represent  the  longitudinal  spectrum 
^^(k)  *  The  Davenport  spectrum  (Ref's.  9  or  10)  is  used  extensively  in  the  building 

aerodynamics  field,  but  it  predicts  that  at  low  frequencies,  <J>uu(k)oCk.  This  is 
inconsistent  with  the  notion  that  <t>^(k)  approaches  a  constant  value  proportional 
to  the  integral  scale  at  low  frequencies,  and  this  model  thus  departs  from  the 
data  at  low  frequencies.  A  modification  of  Davenport's  spectrum  which  removes 
this  objection  has  been  suggested  by  Harris  (Ref's.  58  or  62)  but  this  turns  out 
to  be  identical  to  the  von  K£rm£n  spectrum  below.  The  isotropic  Dryden  spectrum 
corresponding  to  Eq.  (55)  for  the  vertical  component  is  given  by 

h  u2  1  * 


u 


(27TLuxk)‘ 


(60) 


but  like  its  counterpart,  it  is  rejected  due  to  the  -2  slope  it  predicts  for  large 
k.  The  von  K£rm£n  model  is 


<t>  (k)  = 

uu'  ' 


[1  +  70.7(Luxk)2]  5/6 


(61) 


and  it  of  course  does  not  suffer  from  this  objection.  This  is  also  the  case  for 
the  model  spectrum  suggested  by  Fichtl  and  McVehil  (Ref  .39),  similar  to  that 
suggested  for  the  lateral  component,  which  fits  the  data  they  obtained  for  neutral 
and  unstable  conditions  up  to  z  ~  500  ft.  This  equation  is,  non- dimensionally, 


k  <t>  (k)  C  f/f  U 

uu  u  '  p 

o  "  2  u  ru  5  r  / 3  (62) 

pu  u;  [1  +  1.5(f/fpU)  ] 

where  Cu  and  r^  depend  on  stability  and  Pu  and  fpU  depend  on  height  and  stability. 

Finally,  Berman  (Ref.Uo)  has  presented  graphs  which  can  be  used  to  obtain  spectral 
values  for  any  frequency  as  a  function  of  height  and  for  different  stabilities, 
for  heights  up  to  z  -  4 50  ft. 

The  results  of  Gunter  et  al  (Ref.^3)  above  the  surface  layer  show  that, 
as  for  the  lateral  and  vertical  component  spectra,  the  von  K£rm£n  model  provides 
an  excellent  fit  to  the  data  for  various  stabilities  and  surface  conditions,  if 
a  suitable  scale  value  is  used  in  the  model.  The  curves  of  Fig.  6  show  that  it 
also  fits  the  data  well  in  the  surface  layer,  at  least  for  larger  heights.  The 
discrepancy  between  the  results  of  Elderkin  at  z  -  10-20  ft.  and  the  other  two 
curves  emphasizes  the  difficulty  in  modelling  the  spectra  due  to  the  lack  of 
dynamic  similarity. 

Figure  6  also  shows  horizontal  bars  to  indicate  the  range  of  natural 
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frequencies  inrpurtant  to  aircraft  and  structures.  These  bars  have  been  adapted 
from  those  of  lenrplin  (Ref.30)  and  are  plotted  for  a  representative  height  of 
z  *  200  ft.  For  the  'tall  buildings  and  bridges'  bar,  a  mean  velocity  of  50  fps. 
has  been  assumed  for  a  typical  fnsquency  range  of  from  0.1  to  1  cps. 

It  is  to  be  noted  that  none  of  the  model  spectra  for  the  three  velocity 
components  takes  account  of  the  viscous  subrange  in  which  ^^(kJoC  k“?  .  This  is 

in  general  not  a  serious  omission  since  this  region  is  usually  beyond  the  frequency 
range  that  affects  the  motions  of  aircraft  or  the  oscillations  of  exposed  structures. 
However,  as  Templin  shows  in  Ref.30,  if  simulation  of  the  planetary  boundary  layer 
is  attempted  on  too  small  a  scale,  the  viscous  subrange  may  well  be  near  the  region 
of  interest,  resulting  in  significant  errors  if  it  is  ignored.  This  point  should 
thus  be  considered  when  the  scale  of  an  experiment  is  being  determined. 

3.6.4  Cross-Spectra 

In  general,  cross-spectral  data  is  available  only  for  the  different  velo¬ 
city  components  at  a  single  point,  but  as  for  the  power  spectral  densities,  it  is 
only  these  spectra  that  are  of  major  interest.  For  the  $uv(k)  and  $  (k)  spectra, 

Gunter  et  al  have  obtained  measurements  at  z  =  250  ft.  and  750  ft.  for  various 
roughness  and  stability  conditions  and  have  presented  the  results  in  the  form  of 
the  coherences  7uv^(k)  and  7yw^(k).  As  pointed  out  in  Sec.  3 -3 - 3 >  these  coherences 

are  generally  less  than  ~  0.15,  and  tend  to  be  smallest  for  k  >  ~  2/z.  Elderkin 
(Ref. 17)  has  measured  the  co-spectral  parts  of  the  cross-spectra,  Cuy(k)  and  Cvw(k) 

for  heights  up  to  260  ft.  in  various  stabilities.  His  results  are  given  typically 
by  Fig's.  7  and  8  and  while  these  co-spectra  are  not  in  general  significantly 
smaller  than  the  Cuw(k)  co-spectra,  the  values  tend  to  be  centered  about  the  zero 

axis  in  most  stability  cases.  Also,  these  co-spectra  do  not  become  truly  zero  for 
k  less  than  about  3/z,  with  the  values  becoming  larger  for  lower  values  of  k.  As 
for  the  quadrature  component  spectra  Q  (k)  and  G(k),  these  are  generally  assumed 
to  be  zero. 

The  <t>  (k)  cross-spectrum  differs  considerably  from  the  other  two  cross¬ 

spectra  above.  The  quadrature  component  Quw(k)  is  generally  quite  small,  but  its 

sign  is  useful  for  indicating  what  part  of  a  moving  eddy  with  transverse  vorticity 
contributed  to  the  reading.  Thus,  very  qualitatively,  the  height  of  the  centre 
of  gravity  of  an  eddy  having  a  particular  frequency  k  can  be  estimated  by  knowing 
the  sign  of  Quw(^0)  at  a  few  values  of  the  height  (se8  Ref .4),  As  for  the  co¬ 
spectrum  C  (k),  it  is  negative  at  low  frequencies,  with  the  absolute  value  de¬ 
creasing  toward  zero  for  increasing  frequency.  Elderkin' s  data  show  that  it 
reaches  zero  at  about  the  same  frequency  as  do  the  Cuy(k)  and  C^w(k)  co-spectra, 

and  thus  while  the  small  eddies  may  have  considerable  energy,  they  do  not  contribute 
to  the  correlations  between  the  velocity  components.  This  result  has  been  known 
for  some  time,  and  is  consistent  with  the  concept  of  local  isotropy.  Elderkin 
also  finds  that  as  for  the  power  spectra,  there  is  a  slight  shift  of  cuw(k) 

ward  higher  frequency  with  increased  stability.  In  addition,  Elderkin' s  plot  of 
this  co-spectrum  in  similarity  coordinates  shows  that,  like  the  vertical  velocity 
power  spectrum,  it  displays  fairly  reasonable  dynamic  similarity.  Indeed  the  cross 
spectrum  $  (k)  should  not  vary  too  much  with  altitude  through  the  surface  layer 

over  homogeneous  terrain,  since  the  area  under  it  is  uw  which  is  assumed  constant. 


Analytical  representation  of  the  (k)  cross  spectrum  has  been  discussed 

uw 

by  Lumley  and  Panofsky  (Ref .4)  on  the  bt.sis  of  the  similarity  of  the  co-spectrum 
Cuw(k) .  They  have  suggested  an  expression  for  Cuw(k),  but  unfortunately  this  does 

not  fit  the  data  at  very  low  frequencies  or  in  the  inertial  subrange.  A  useful 
model  for  $uw(k)  has  been  suggested,  however,  by  Case  et  al  (Ref.  53),  in  which 
Q  (k)  is  assumed  zero  such  that  <t>  (k)  =  C  (k).  They  used  data  taken  near  the 

Uw  UVf 

wal?.  in  a  wind  tunnel  boundary  layer  and  also  in  the  surface  boundary  layer  over 
the  ocean  to  obtain  the  root-coherence  functions  defined  from  Eq.  (28)  by 


It  was 


y  (k)  ~ 

uw'  ' 


|  4>  (k)  1 

'  uw'  '  1 _ 

[<t>  (k)d>  (k)]1/2  * 

L  uu  ww  J 


;  ound  that  the  expression 


7  (k)  = 

uw'  ' 


[1 


0.395(L„xk)2]1/2 


u 


(63) 


(64) 


provided  a  good  representation  of  the  experimental  data.  The  quantity  7q  is  a 
constant  whose  value  is  determined  by  the  power  spectra  <^>uu(k-)  and  <t>ww(k)  and 
the  value  of  the  Reynolds  Stress  coefficient  uw/u'w' .  That  is,  from  Eq.  (24), 


<t>  (k)  dk 

uw'  ' 


$  (k)  dk 

uw'  ' 

since  if  Q  (k)  is  zero,  $  (k)  =  C  (k)  =  C  (-k).  Also,  since  C  (k)  is  nega- 

uw  uw  uw  uw  *  uw  '  ° 

tive  at  all  values  of  k  until  it  reaches  zero,  we  may  write 


uw  *  -  /  l<*>uw(k)  I  dk* 


Thus  from  Eq.(63),  we  get 


_  r 00  1/2 

uw  =  -  I  7  (k)[  <t>  (k)$  (k)l  dk 

uw'  ' 1  uu  '  ww  ' J 


and  using  Eq.  (64), 


uw  -  -7 


•  00  r  $  (k)  <t>  (k) 

uu  ww 1 

O  1  +  0.395(L  xk)‘ 


1/2 


dk. 


Now  recall  from  Eq.  (r>?)  that  in  the  surface  layer, 


(65) 


30 


—  Ip? 

UW  =  -  0.31  U'  w'  -  -0.31  V u  *w  , 

such  that  when  this  is  combined  v;it.h  Ec.  (65) ,  we  get 
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(66) 


dk 


Thus  7  can  be  found  for  any  choice  of  the  <t>  (k)  and  <t>  (k)  spectra,  and  the 

'o  uu  ww 

cross  spectrum  is  obtained  from  Eq’s.  (63)  and  (64)  as 


r  $  (kVt>  (k) 
.  \  i  uu'  ww 

<J>  (k)  =  -7 


UW 


O  ! 


L  1  +  0.395(L„xk)2 
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(67) 


Evaluation  of  7  as  a  function  of  2  has  been  done  using  the  von  K^rm^n  spectra 
(Sec. 1^.5  and  Fig. 12). 

3.6.5  Cross-Correlation  Spectrum 

Consider  as  a  special  case  of  the  general  correlation  tensor  of  Eq.(3) 


the  tensor 


R..(r  ,r  )  =  R.  .(r  ,r  ,0,0). 
zy  x®  y®  x*  y * 


(68) 


In  order  to  obtain  a  tensor  of  true  two-dimensional  spectrum  functions,  one  would 
Fourier  transform  this  tensor  twice,  such  that 


.(k  ,k  )  «  4  1 1  R. ,(r  ,r  ) 
ij  x*  y'  JJ  ±y  x*  yJ 


-i£7r(k  r  +  k  r  ) 
x  x  y  y 

e  dr  dr  . 

x  y 


(69) 


Similarly,  if  the  correlations  R.  .(r  )  are  Fourier  transformed  once,  one  gets  a 

1J  X 

tensor  of  true  one- dimensional  spectrum  functions 


0.  .(k  )  =  2  /  R.  .(r  )  e 

r  n  n  '  v '  1  in'  v' 


-1277k  r 

X  X 


1J  X 


X' 


dr 


as  in  Eq.  (20a).  However,  suppose  now  that  the  two-dimensional  correlation  functions 

R.  .(r  ,  r  )  are  transformed  only  cnce,  with  respect  to  r  .  In  this  case,  the 
x*  y'  j  ^  x 

results  are  not  true  spectrum  functions  in  the  sense  defined  above,  but  neither 

are  they  correlation  functions.  The  result  of  such  a  transformation  is 
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Assuming  Taylor's  Hypothesis,  k  may  be  replaced  by  k  -  n/U,  and  since  r  =  Ut , 

X  X 
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(k,  r  )  =  2  0  /  <R.  .(r  ,t) 

ij  y  J_ y 
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(71) 


for  homogeneous  turbulence.  Indeed,  if  the  turbulence  is  not  homogeneous  in  the 
y-direction,  this  quantity  is  a  function  of  the  two  points  y  and  y^  which  are 
separated  by  r  rather  than,  of  only  r  .  That  is, 


r -i277nT 

^(k,  y]3  y?)  =  2  uj  ^ij(y1>  y2>T)  e 


This  quantity  can  also  be  expressed  in  terms  of  a  coherence  function,  defined  by 
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where  (k,y^)  refers  simply  to  the  power  spectrum  tfh^k)  of  the  i-component 

obtained  at  y.  ,  and  similarly  for  «t>^(k,  y^).  Physically,  ^^(k,  y^,  y2)  may  be 

interpreted  as  the  zero-time-delay  cross-correlation  (or  mean  product)  of  the 
parts  o'"  the  two  velocity  components  i  and  j  contained  in  their  respective 
spectra  at  the  frequency  k,  and  this  has  led.  Davenport  (Refs. 9  and  10)  to  use 
the  term  'narrow-band  cross-correlation',  of  which  7„(k,  y^,  y2)  as  defined 

here  is  the  magnitude.  It  may  also  be  interpreted  as  the  frequency  distribution 
of  the  spatial  cross- correlation  of  the  i-  and  j-components,  leading  to  the  term 
'cross-correlation  spectrum'  also  used  by  Davenport.  Finally,  it  is  pointed  out 
that  since  .(k,  y.  ,  y  )  is  a  form  of  cross-spectrum,  it  may  be  expressed  as 

1 J  X  f— 

the  r.uia  of  real  and  imaginary  parts.  That  is, 

Vk’  yi’  ya>  *  c«<k>  V  y2>  T  1  Vk>  yi>  ys> 

and  thus  Eq.  (72)  becomes 
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In  the  case  of  homogeneous  turbulence,  =0  and  the  y^  and  y^  dependence  is 
replaced  by  such  that 
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Davenport  measured  the  quantity  7uu(k,  r  )  as  defined  above  for  the 

longitudinal  velocity  components  at  two  points  separated  laterally  across  a  strong 
turbulent,  flow  in  a  wind  tunnel  (see  Ref's.  9  or  22).  His  results  showed  that 
for  various  values  of  r  ,  the  data  collapsed  on  a  curve  given  by 
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where  c  ^  8,  For  the  turbulence  in  the  atmosphere,  measurements  have  shown  that 
for  separations  in  the  vertical  direction,  cross  correlations  of  the  longitudinal 
component  indicate  the  existence  of  a  similar  relationship.  That  is,  if  the  sep¬ 
aration  in  the  y-di.rection  as  above  is  replaced  by  separation  in  the  z-direction 
and  the  coherence  is  defined  by 


Cl/  ^ 


(75) 


then  measurements  at  various  vertical  spacings  have  shown  that 
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where  the  coefficients  c1Q  and  c^  are  defined  depending  on  whether  the  velocity  0 
at  height  z  or  at  z  =  10  m.  is  used.  Generally,  the  value  0^  is  preferable 
here,  since  the  increase  of  U  with  height  means  that  U  is  different  at  z^  and  z0 

and  thus  the  question  arises  as  to  which  value  of  0  to  use.  Davenport  in  Ref,  10 
quotes  values  of  c^Q  ranging  from  5  to  10  for  heights  from  300  to  830  ft.  over 
various  types  of  terrain.  The  coefficient  appears  to  be  independent  of  height, 
but  seems  to  be  generally  smaller  for  rough  surface  conditions  such  as  urban  areas 
than  for  smooth  surfaces. 


Surry  (Ref  .12)  investigated  grid  turbulence  in  a  wind  tunnel  and  found 
that  a  modified  version  of  the  Davenport  relation  of  Eq,  (74)  better  fits  his 
data  for  7uu(k>  ry)  than  the  original  expression.  This  modified  empirical  relation 

-(ckr  )a 

7  (k,  r  )  -■  e  (77) 

uu  *  y' 

where  c  «  6.4  and  a  *  1.4.  In  addition,  Surry  points  out  that  a  curve  of  this 
form  with  a  =  1.4  and  c  =  8  provides  a  better  fit  to  Davenport's  own  data  than 
Eq.  (74).  Finally,  Harris  (Ref. 58)  has  shown  that  the  data  is  better  represented 
by  a  more  complex  analytical  expression  based  on  homogeneous  isotropic  theory. 

It  can  be  shown  for  the  von  K£rm£n  model  that 

5/6  11/6 
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33 


where  ti  =  2rr  U  „  \/2  +  It  /LU  and  n  =  nL/0,„.  The  parameter  L  is  a  length  scale 

y  10  Y  '  '10  _ 

and  K  ,,  and  K  ,r  are  modified  Bessel  functions  of  the  second  kind.  If  2, 

5/6  1/6 

which  is  the  case  for  most  of  the  frequency  range,  the  above  expression  for  tj 

becomes 

T)  ~  27mr y/U  =  ^7rkry 


and  direct  comparison  of  Eci.  (78)  with  the  above  empirical  relations  becomes  quite 
simple.  This  comparison  is  shown  in  Fig.  9>  along  with  experimental  results  meas¬ 
ured  in  wind  tunnel  grid  turbulence  by  Bearman  and  quoted  by  Harris  in  Ref.  58. 

It  is  seen  that  the  theoretical  expression  is  superior  since  it  accounts  for  the 
negative  portions  of  the  curve,  while  the  empirical  curves  do  not.  Indeed,  neither 
Surry's  rcr  Davenport's  experimental  data  showed  negative  results,  presumably 
because  they  found  results  only  for  T)  up  to  ~  2.2. 


As  for  atmospheric  data,  Harris  compares  data  measured  at  points  separated 
vertically  above  the  surface  with  the  corresponding  theoretical  expression 

7uu(k>  rz>  =  {  (  2  )  K5/6^>  -  (  2  )  K1/6(T,)  } 

where  n  =  27rkr  .  His  results  show  very  good  agreement  for  heights  above  the  surface 
z 

layer,  but  deterioration  of  this  agreement  within  the  surface  layer.  This  result 
is  not  at  all  surprising,  since  it  is  an  indication  of  the  closer  adherence  of 
atmospheric  turbulence  to  isotropic  relations  above  the  surface  layer  than  in  it, 
and  this  feature  is  indicated  by  virtually  all  the  other  turbulence  data. 

3.7  Integral  Scales  of  Turbulence 

Integral  scales  of  turbulence  have  been  defined  in  Sec.  2.2  and  of  the 
nine  scales  defined  by  Eqs.(l2),  the  scales  L  x  and  I*wx  of  Eq.  (12(a))  are  of 

the  most  interest.  In  order  to  obtain  estimates  of  these  scales  from  experimental 
data,  three  basic  methods  can  be  used.  The  first  of  these  is  what  might  be  called 
the  direct  method,  since  it  involves  obtaining  the  actual  correlation  curve  in 
the  scale  definition  and  integrating  to  obtain  the  area  under  it.  In  the  case 
of  the  above  scales,  the  single  point  time  correlation  curve  is  generally  measured 
and  Taylor's  Hypothesis  assumed  to  convert  the  time  scales  so  obtained  to 
length  scaler.  Thus,  for  example,  from  Eq.  (l6) 
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The  difficulty  that  often  occurs  with  this  approach  is  that  at  the  low  frequency 
(large  t  or  spatial  separation)  end,  the  correlation  does  not  approach  zero 
asymptotically.  Instead,  due  to  the  existence  of  very  slow  variations,  the 
correlation  obtained  from  a  finite  recording  period  has  undulations  about  the 
zero  axis  at  the  low  frequency  end  (viz  .Ref's.  24  and  48).  One  attempt  at  avoid¬ 
ing  this  problem  is  the  use  of  the  second  method  for  obtaining  scales.  In 
this  method,  one  uses  for  the  scale  the  value  of  the  space  (or  time)  separation 
at  which  the  normalized  correlation  curve  has  dropped  to  a  value  of  l/e.  Thus 
assuming  the  correlation  is  a  true  exponential,  such  that,  for  example, 
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as  in  Eq.  (5*0 »  then  the  above  value  of  is  the  true  scale  .  As  one  would 

expect,  however,  this  ic  usually  not  the  case,  and  some  form  of  correction  to  this 
value  is  generally  required,  This  correction  usually  requires  the  application  of 
the  first  method  at  least  a  few  times. 


The  third  method  for  obtaining  integral  scales  is  from  the  one-dimensional 
power  spectra.  In  all  the  model  spectra,  there  is  a  free  scale  parameter  which 
locates  the  curve  on  the  frequency  axis,  and  for  the  von  K^rmfm  or  Dryden  models, 
this  parameter  is  the  true  integral  scale  L  x.  For  any  particular  model,  this 
scale  parameter  can  be  related  to  the  value11  k  *  at  which  the  non-dimens ionalized 

spectrum  has  its  maximum  value.  For  example,  From  the  von  K^rmfin  model  for  the 
longitudinal  conponent  spectrum,  Eq,  (6l),  we  get 

Lux  -  0.1*46/kpu  (79) 

while  the  Dryden  model  of  Eq.  (60)  gives 

LuX  -  l/27rkpu  =  0.159ApU.  (80) 


For  the  vertical  component  models,  Eq.  (56)  yields 


and  from  Eq.  (55) 


LUX  -  0.212/kp" 


LUX  -  0.23V*p". 


(01) 

(82) 


The  lateral  conponent  models  will  of  course  give  the  same  expressions  as  Eq's. 
(8l)  and  (82)  with  kpW  replaced  by  kpv.  Suppose  now  that  either  the  von  KkrmCn 

or  the  Dryden  mouel  is  assumed  to  fit  the  experimental  spectral  data  for  a 
particular  component.  The  scale  can  thus  be  simply  obtained  from  the  location 
k  *  of  the  peak  of  the  experimental  curve  by  using  the  appropriate  expression 
aSove.  This  method  for  obtaining  scale  estimates  is  used  quite  commonly,  but 
it  has  two  major  drawbacks.  In  the  first  place,  the  peak  of  the  experimental 
data  is  not  always  clearly  defined,  particularly  for  the  longitudinal  and  lateral 
components,  and  it  can  therefore  be  difficult  to  obtain  a  good  estimate  of  k  *. 

In  addition,  the  peak  is  generally  located  near  the  low  frequency  end  of  the 
spectrum,  where  the  data  is  not  as  reliable  as  it  is  at  larger  frequencies,  and 
further  errors  may  be  introduced. 


In  obtaining  scale  values  by  the  third  method  outlined  above,  con¬ 
siderable  care  must  be  exercised  in  interpreting  exactly  what  has  been  obtained. 
As  previously  stated,  both  the  von  Karmen  and  the  Dryden  spectral  models  are 
for  isotropic  turbulence,  and  the  vertical  and  lateral  component  models  are 
written  in  terms  of  the  longitudinal  component  scale  Lux  on  the  assumption  that 

the  isotropic  relation  L  x  =  2L  x  =  2L  x  is  valid.  However,  if  the  turbulence 
is  not  isotropic,  the  isotropic  models  may  still  fit  the  experimental  data  quite 
well,  but  the  above  scale  relation  is  not  valid.  In  particular,  this  is  the  case 
in  the  lower  atmosphere.  Thus  if  one  applies  Eq,  (8l)  or  (82)  to  experimental 
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data  for  the  vertical  velocity  component,  it  is  in  fact  not  that  has  been  ob¬ 
tained,  but  rather  2  L  X,  and  from  this  data  one  can  obtain  information  only  about 

the  vertical  component,  not  the  longitudinal.  This  point  explains  why  Taylor 
^Ref .51) ,  for  example,  concludes  that  Lux  is  proportional  to  height  z  from  the 
surface  up  to  z  =  1C00  ft.  when  moot  results  indicate  a  much  weaker  variation 
with  height.  Taylor's  results  are  based  on  vertical  component  data  taken  from 
Saunders(Ref .20)  and  Panofsky  and  McCormack  (Ref. 29)  and  fitted  with  the  von 
K^rmfc  vertical  component  model.  The  data  shows  that  kpw  is  proportional  to 

height,  such  that  while  Eq.  (8l)  indicates  L^X  should  thus  be  proportional  to 
height,  it  is  in  fact  L  x  that  is  proportional  to  height.  No  conclusions  can  be 

x  W  XX 

made  about  L  except  perhaps  at  larger  heights  where  =  2Lw  is  a  more  reason¬ 
able  approximation.  Probably  the  best  way  to  avoid  confusion  here  when  using  the 
von  KarmCo  or  Dryden  models  is  to  replace  L^x  by  2L^X  in  the  lateral  component 
equation  and  by  2L  x  in  the  vertical  component  equation.  Thus  for  the  von  Karman 
model,  Eq's.  (79)  and  (8l)  become 

L  X  =  0.lU6/kpU,  LwX  =  0.106/kpW  (83) 

and  for  the  Dryden  model,  Eq's.  (80)  and  (82)  become 

Lux  =  0.159/kpu,  Lwx  =  0.117/kpw.  (84) 

x 

In  the  literature,  the  bulk  of  the  data  available  is  for  the  scales 
and  LwX.  As  one  might  expect,  these  scales  are  functions  of  height,  thermal 

stability,  and  surface  conditions  since  physically  they  represent  the  size  of  the 
predominant  eddies  in  the  turbulence.  As  with  the  longitudinal  and  lateral  com¬ 
ponent  variances,  the  effect  of  surface  conditions  is  not  adequately  expressed  by 
roughness  length,  since  large  scale  non-uniformities  such  as  hills  and  mountains 
have  a  significant  influence  on  the  scales.  For  example,  the  results  of  Berman 
(Ref.40)  shew  little  effect  of  roughness  length  on  scale  for  uniform  terrain 
while  Lappe  (Ref .28)  indicates  a  distinct  increase  in  scale  with  increasing  sur¬ 
face  non- uniformity.  Lappe  also  shows  that  the  effect  of  the  surface  is  felt 
most  strongly  at  lower  heights,  end  that  as  height  increases  toward  1000  ft., 
terrain  type  becomes  less  important.  This  tendency  was  also  found  by  Gunter  et 
al.  (Ref. 43)  whose  results  show  that  increasing  surface  non-uniformity .increases 
scale  at  z  =  250  ft.  but  has  little  effect  at  750  ft.  Conversely,  an  increase 
of  scales  with  decreasing  stability  was  found  by  Gunter  ot  al  to  be  quite  noticeable 
at  z  =  750  ft.  but  not  at  250  ft.  To  summarize,  then,  it  may  be  concluded  that 
increasing  surface  non-uniformity  tends  in  general  to  increase  integral  scales, 
but  more  so  at  small  heights  than  at  large;  decreasing  thermal  stability  tends  to 
increase  scale,  but  more  so  at  large  heights  than  at  small.  Also,  roughness 
length  has  little  effect  on  scales,  while  increasing  large  scale  non- uniformities 
of  the  surface  tends  to  increase  them. 

As  for  the  effect  of  height  on  scale,  it  is  generally  agreed  that  all 

of  L  X,  L  x  and  L  x  increase  with  height,  although  not  necessarily  in  the  same 
u  ’  v  w 

manner.  The  variation  of  the  scale  Lwx  of  the  vertical  component  is  the  one 
which  is  most  consistently  agreed  upon.  Figure  10  shows  the  results  of  most  of 
the  investigations  of  recent  years,  for  neutral  air  over  flat  terrain  unless 
otherwise  noted.  In  all  cases,  the  scales  were  obtained  by  the  spectral  fit 
method,  and  in  the  few  cases  when  only  the  spectra  were  given,  Eq.  (83)  was  used 
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to  determine  L  *.  It  is  noted  also  that  for  the  data  of  Ref's.  20,  21  and  28,  the 
published  values  for  L  '  were  divided  by  2  and  used  as  Lwx  since  these  values  were 
determined  from  vertical  component  data.  It  is  evident  that  the  scale  LwX  is  in 
fact  proportional  to  height  as  suggested  previously,  and  a  value  of  0,4o  for  the 
constant  of  proportionality  places  the  curve  in  the  centre  of  the  data  within  a 
reasonable  amount  of  scatter.  Thus  it  is  suggested  that  the  relation 

LWX  =O.I(z.  (85) 

gives  a  reasonable  estimate  of  this  scale  for  heights  up  to  1000  ft.  in  neutral  air 
over  relatively  flat  terrain.  It  is  also  pointed  out  that  if  one  considers  only 
the  region  outside  the  surface  layer  in  Fig.  10  (i.e.,  fc  =  200  or  300  ft.  to 
z  =  1000  ft.),  the  variation  of  L  x  with  height  appears  to  weaken  somewhat,  as 
suggested  by  Busch  and  Panofsky  (ftef.  hi).  The  additional  line  shown  in  this  region 
is  given  roughly  by 

X  0.73 

Lw  *  2.1  z  (86) 

and  this  may  fit  the  data  slightly  better  in  this  region  than  the  L  X  *  0.4  z 
relation. 


x 

The  results  available  to  date  for  the  scale  Lu  at  various  heights  in 

neutral  air  are  shown  in  Fig.  11.  Here  again,  most  of  the  data  were  obtained  by 
fitting  spectral  models  to  experimental  data  for  the  longitudinal  component.  It 
is  obvious  that  there  is  considerably  greater  scatter  in  this  data  than  in  those 
for  L  x.  However,  some  definite  conclusions  can  be  drawn.  First,  it  can  be  seen 
that  ¥here  is  a  distinct  increase  in  L  x  with  height  over  the  entire  range  0-1000 
ft.  Also,  this  increase  with  height  is  weaker  than  that  for  Lwx,  with  Berman 

(Kef  .40)  suggesting  Lux  oc  z°*2^and  Webb  concluding  LuXoC  z0*^.  Thirdly,  the  data 

for  the  surface  layer  is  considerably  more  conclusive  than  that  for  the  remainder 
of  the  planetary  layer.  In  the  0-200  ft.  range,  most  of  the  data  falls  within 
the  shaded  band  shown.  The  centre  of  this  band  is  given  by  the  line 

X  0*47 

Lux  =  22.7  z  (87) 

and  the  extremities  of  the  band  axe  roughly  50-6O)E>  away  from  this  line.  Since 
this  relation  suggests  more  accuracy  than  is  reasonable  considering  the  scatter, 
the  relation 

LuX  =  20  \TT  (88) 


is  also  shown  in  Fig.  11,  and  it  is  seen  to  reasonably  indicate  the  location  of 
the  band.  Thus  it  is  suggested  that  Eq.  (88)  be  used  to  estimate  L  x  in  the  sur¬ 
face  layer  under  these  conditions,  realizing  that  variations  of  the  size  shown 
are  unavoidable  due  to  the  complexity  of  the  atmosphere.  Above  200  ft.,  Fig.  11 
shows  that  Eq.  (88)  can  also  be  used  with  a  scatter  in  the  data  similar  to  that 
below  200  ft.  However,  it  is  to  be  noted  that  over  the  200-1000  ft.  range,  a 
stronger  or  weaker  variation  could  be  the  case,  as  the  data  are  not  too  conclusive. 
In  particular,  the  line  L  x  ~  0.8  z  is  seen  to  fall  within  the  scatter  band,  as 


is  the  line 


-  4 


0.73 
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(89) 
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Both  these  lines  represent  the  isotropic  relation  Lu  =■  2Lw  based  on  the 
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relations  of  Eq’s.  (85)  and  (86)  for  L  X  in  the  200-1000  ft.  range.  However,  the 

X  ^ 

line  L  =  0.8  z  is  seen  to  be  completely  unacceptable  below  z  =  200  ft.  Thus, 
here  again,  as  from  other  observations,  it  is  seen  that  if  any  form  of  isotropy 
exists  in  the  lower  atmosphere,  it  must  certainly  be  outside  the  surface  layer. 

In  Fig.  11,  the  fact  that  LuX  increases  with  height  is  indicated  by 

most  of  the  data  shown,  including  that  of  Davenport  (Ref .37).  This  is  contrary 
to  the  conclusion  of  Elderkin  (Ref. 17)  and  Fichtl  and  McVehil  (Ref. 39)  who  state 
that  Davenport's  data  showed  the  peak  of  the  longitudinal  spectrum,  and  thus  the 
scale  L  ,  to  be  invariant  with  height.  However,  the  invariance  of  this  data 
with  height  results  only  from  the  fact  that  the  variable  n/U^Q  was  used  as  the 

abcissa,  where  U,  Q  is  the  mean  velocity  at  z  =  10  m.  Indeed,  as  Davenport  himself 

points  out  in  Ref. 10,  if  n/U  were  used  as  the  independent  variable,  where  U  is  the 
mean  velocity  at  height  z,  then  the  location  of  the  peak, which  would  now  be  k^u 

as  previously  defined,  would  indeed  show  a  variation  with  height.  Thus  the  'scale' 
I,  used  by  Davenport  in  Ref .10  (and  which  was  equal  to  1200  m.  in  Ref.  37)  and 
also  by  Harris  (Ref's.  58  and  62)  is  constant  with  height,  but  only  because  it 
was  defined  using  n/U^.  By  comparing  Harris'  spectral  model  with  that  of  von 
Kerman,  this  scale  may  be  found  in  terms  of  the  true  integral  scale  as 


and  it  is  thus  easily  seen  that  since  L  is  constant, 

Luv  o/o10*  *a 

using  Eq.  (44).  Thus  Davenport's  results  do  not  indicate  constant  longitudinal 
scale  but  rather  that  this  scale  increases  with  height  in  the  same  way  as  the 
mean  velocity  does. 


x 

Very  little  data  is  available  for  the  scale  ,  and  few  definite 

conclusions  can  be  drawn.  Gunter  et  ai  (Ref. 43)  found  that  oversAllsStatilities 
and  terrain  types,  the  average  value  of  Lvx  was  180  ft.  at  z  =  250  ft.  and 

230  ft.  at  z  =  750  ft.,  these  values  being  roughly  equivalent  to  those  found 
for  1  x  and  half  the  values  found  for  L  x.  The  spectra  obtained  by  Fichtl  and 
McVehil.  (Ref. 39)  were  used  to  find  thatuin  neutral  air. 


0.42 

35.2  z 


x 

for  z  ~60-500  ft.,  indicating  a  variation  with  height  similar  to  that  for  Lu 
but  yielding  values  that  are  roughly  the  same  in  magnitude  as  L  x  rather  than 

fj  x.  Elderkin1  s  data  for  z  =  10-20  ft.  yield  L  x  z,  such  thkt  L  x  is  larger 
w  v  v 

than  the  values  of  I.  x  obtained  from  his  data  but  smaller  than  the  L  x  values 

w  u 

so  obtained.  Thus  about  all  that  can  be  said  about  L  x  is  that  it  increases 

v  x 

with  height,  and  above  the  surface  layer  it  is  probably  roughly  the  same  as  L  . 


As  for  the  remaining  scales  defined  by  Eqs.(l2),  little  information  is 

available.  Panofsky  (Ref .24)  gives  some  information  based  on  correlation  curves 

measured  over  smooth  terrain  (z  =1  cm)  at  z  =  6.5  m.  From  so-called  'semi-scales' 

o 
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obtained  from  these  correlations,  he  found  that  in  stable  air,  L,,  ~  8  L  y,  a 

11  u 

result  which  was  also  found  by  Grant  (Ref .2$)  in  a  wind  tunnel  boundary  layer. 
Thus  the  eddies  are  elongated  in  the  mean  wind  direction  in  stable  air, and  also 
presumably  in  neutral  air.  The  lateral  component  scale  was  also  found  to  be 
elongated  in  the  mean  wind  direction  in  stable  air,  with  L  x  ~  3  L  Y,  and  also 
L  y  ~  L  y.  In  unstable  air,  however,  Panof sky's  results  Xhow  an  entirely 

V  Vi  JC 

different  situation,  in  that  Lu  was  found  to  be  increased  by  a  factor  of  about 
six  over  the  stable  air  value,  and  all  scales  were  roughly  the  same  size,  or 

Lx  ~lY~Lx~LY.  Finally,  Panofsky  also  found  L  x  and  L  Y  much  smaller 
uuvv  *  J  w  w 

than  the  other  scales  at  this  height  for  both  stable  and  unstable  air,  a  result 

not  surprising  considering  the  nearness  of  the  surface.  Thus  the  anistropy  of 

the  turbulence  very  near  the  surface  is  emphasized  by  these  results,  in  addition 

to  the  strong  effect  of  stability.  Panofsky  sums  up  his  findings  by  stating  that 

on  a  clear  day  (i.e.,  unstable  air)  the  scales  L  x,  L  Y,  L  x,  and  L  Y  are  all 
'  '  ’  '  u  u  v  v 

large  and  of  the  same  order  of  magnitude,  while  at  night  (stable  air)  the  scales 

are  much  smaller  and  that  of  the  u-component  is  much  smaller  across  the  wind  than 

along  it.  For  a  further  qualitative  description  of  the  eddy  structure  in  the 

atmosphere,  the  reader  is  referred  to  Ref.  4,  page  210. 

IV.  MATHEMATICAL  MODEL  OF  THE  PLANETARY  BOUNDARY  LAYER 


As  a  summary  to  the  preceding  description  of  flow  characteristics  in 
the  planetary  layer,  a  mathematical  model  is  suggested,  based  on  certain  simpli- 
fying  assumptions.  This  model  is  intended  to  be  a  compromise  between  precision 
and  simplicity,  and  the  validity  of  any  part  of  it  can  be  determined  by  referring 
to  the  appropriate  section  of  the  preceding  discussion. 

4.1  Assumptions 

(1)  The  mean  wind  velocity  is  assumed  to  be  strong  enough  such  that  the 
atmosphere  may  be  considered  neutrally  stable. 

(2)  Taylor's  Hypothesis,  or  the  assumption  of  frozen  flow,  is  assumed  to 

be  valid  over  the  entire  gust  spectrum  range  of  interest  (A  =  0/n ~ 50-7000  ft). 

(3)  The  gross  features  of  the  terrain  are  assumed  to  be  relatively  uniform, 
such  that  roughness  length  zq  represents  the  surface  effects  adequately. 

(4)  The  flow  is  assumed  stationary  for  periods  at  least  as  long  as  the 
record  length  used  to  obtain  mean  velocity  values. 

(5)  The  flow  is  assumed  to  be  homogeneous  in  any  horizontal  plane. 

4.2  Mean  Velocity  Profile 

The  mean  velocity  through  the  entire  planetary  layer  ( z  «  0  to  ~  2000  ft,) 
is  given  by 

0/0G  =  (z/z/ 

where  z_  and  a  are  given  in  terms  of  surface  roughness  by  Fig.  1  and  the  gradient 

wind  U_  is  determined  from  local  meteorological  data.  For  the  surface  layer  alone 

G 

( z  *  0  to  ~200  ft.)  the  mean  velocity  is  given  by 
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0/0  =  2.5  In  (z/z  ) 

*  o 

where  0  is  the  friction  velocity. 

1 . 3  Reynolds  Stresses 

The  ratios  of  the  velocity  component  variances  in  the  surface  layer  are 

given  by 


u'/v'/v'/0t  =  2. 5/2. 0/1. 3/1 

and  the  turbulence  intensities  are 

u'  1  v'  0.80  ,  w' 

—  -  -  ,  —  =  -  ,  and  — 

0  ln(z/zQ)  U  ln(z/zQ)  0 


0.52 

ln(z/zQ) 


Above  the  surface  layer,  the  intensities  all  tend  toward  the  same  value,  this  value 
being  determined  by  the  turbulence  intensity  in  the  free  atmosphere  which  is 
assumed  to  be  zero. 


The  Reynolds  Stresses  uv  and  vw  are  assumed  to  be  zero  throughout  the 
planetary  layer.  In  the  surface  layer,  the  uw  Reynolds  Stress  is  determined  by 
the  variance  ratios  specified  above,  such  that  (Eq.52) 

uw  /  u'w'  =  -0.31. 


Above  the  surface  layer,  uw  decreases  in  magnitude  toward  zero. 


4.4  Power  Spectra 


The  von  K^rmfin  equations  are  suggested  as  the  model  for  the  velocity 
component  power  spectra.  The  appropriate  (anistropic)  values  of  integral  scale 
and  variance  must  be  used,  however,  as  specified.  Thus  the  power  spectra  are 
given  by 


uu 


00 


[1  +  70.7(LuXk)2 


(90) 
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’/V 
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1  +  l88.4(2Lyxk)2 

[1  +  70.7(2Lv*k)2  ]  11//6 
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cl'  (k)  =  wL.  hh 
ww  w 


1  +  li8.4(2LwXk)2 
[1  +  70.7(2Lwxk)2]11^ 


(91) 
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(92) 


4.5  Crosa-Spcctra 


The  cross-spectra  (k)  and  <t>  (k)  are  assumed  to  be  zero  in  this  model, 

uv  vw  ' 

along  with  the  Quad- spectrum  ^w( k ) .  The  cross-spectrum  $uw(k)  is  thus  equal  to 
Cuw(k)  and  is  given  in  the  surface  layer  by 


^uw  (k) 


4>  (k)  (k) 

uu  wwv  ' 

1  +  0.395(Luxk)? 


\/? 


(93) 


or  in  non-dimensional  terms, 


7  00  - 

uw'  ' 


[1  +  0~395(Luxk)2  ]1/f" 


(94) 


The  value  of  7q  can  be  determined  from  Eq.  (66)  using  the  power  spectra  of  Eqs .  (90) 
and  (92)  and  the  scale  values  given  in  the  following  section.  The  variation  of 
7Q  with  height  is  determined  completely  by  the  ratio  Lwx/Lux,  and  this  variation  is 

shown  in  Fig. 12  for  LwX/LuX  as  found  from  Eqs. (95)  and  (96),  In  the  figure,  y q  is 
shown  to  be  constant  with  height  above  the  surface  layer,  since  Lwx/Lux  *  0.5  *  a 
constant.  In  fact,  however,  yQ  will  decrease  as  uw/u'w’  decreases,  since,  as  indi¬ 
cated  by  Eq.  (66),  is  proportional  to  this  value. 

4.6  Integral  Scales 


For  the  planetary  layer 
and  Lw  are  given  by 


as  a  whole,  assume 


L  x  *  L  x  and  the  scales  L  x 
v  w  u 


L  X  =  20  •J~z 
u 


(95) 


and 

L  X  =  0.4  z  . 
w 


When  the  region  above  the  surface  layer  is  being  considered  alone,  the  expressions 

L  x  =  4.2  z0,73  (96) 

u 

and 


L  X  =  2  1  z°'73 
w 


are  suggested,  and  thus  in  this  region,  the  isotropic  relation 


LX  =  2LX  =  2LX 
U  V  w 


is  valid. 
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FIG.  1  PARAMETERS  OF  THE  POWER  LAW 

FOR  MEAN  WIND  VELOCITY  (FROM  REF.  9) 


GRADIENT  HEIGHT  Zc  -  FT. 
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FIG.  3  TURBULENCE  INTENSITY  VARIATION  WITH  HEIGHT 
IN  NEUTRAL  STABILITY  OVER  FLAT  TERRAIN 

o: 


FIG.  4  VERTICAL  COMPONENT  POWER  SPECTRA  IN  THE  SURFACE  LAYER 


LONGITUDINAL  COMPONENT  POWER  SPECTRA  IN  THE  SURFACE  LAYER 
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FIG.  7  COSPECTRA  BETWEEN  LONGITUDINAL  AND 
LATERAL  VELOCITY- NEUTRAL  STABILITY 


FIG.  8  COSPECTRA  BETWEEN  LATERAL  AND 
VERTICAL  VELOCITY-STABLE  AIR 
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a  EXPERIMENTAL  POINTS  (BEARM AN  FROM  REF.  58) 


FIG.  9  CROSS-CORRELATION  SPECTRA  FOR  GRID-GENERATED  TURBULENCE 


FIG.  10  VERTICAL  COMPONENT  SCALE  VARIATION  WITH 

HEIGHT  IN  NEUTRAL  STABILITY  OVER  FLAT  TERRAIN 


FIG.  11  LONGITUDINAL  COMPONENT  SCALE  VARIATION  WITH 
HEIGHT  IN  NEUTRAL  STABILITY  OVER  FLAT  TERRAIN 


